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Abstract 

In this paper we consider the Grioli-Koiter-Mindlin-Toupin linear isotropic indeterminate couple stress 
model. Our main aim is to show that, up to now, the boundary conditions have not been completely un¬ 
derstood for this model. As it turns out, and to our own surprise, restricting the well known boundary 
conditions stemming from the strain gradient or second gradient models to the particular case of the inde¬ 
terminate couple stress model, does not always reduce to the Grioli-Koiter-Mindlin-Toupin set of accepted 
boundary conditions. We present, therefore, a proof of the fact that when specific “mixed” kinematical and 
traction boundary conditions are assigned on the boundary, no “a priori” equivalence can be established 
between Mindlin’s and our approach. 
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1 Introduction 


Higher gradient elasticity models are nowadays increasingly used to describe mechanical systems with underlying 
micro- or nano-structures (see e.g. among many others [7J [5UJ 0J EH EH HH 1H1221EH ESI ESI ESI EH 0H EO) or to 
regularize certain ill-posed problems with higher gradient contributions (see e.g. [55; 62, 1221 21-21-155] I. Such 
higher gradient models, together with the more general class of micromorphic models [511 El [Ml EB EH HU EH 0 , 
have been also proved to be a useful tool for the description of micro-structured materials showing exotic 
behaviours in the dynamic regime (see e.g. O [531 EH Ell EH EES Ell 1221 IBS] ) • 

One among such higher order models which was introduced at the very beginning is the so called inde¬ 
terminate couple stress model in which the higher gradient contributions only enter through gradients on the 
continuum rotation. We place ourselves in the context of the linear elastic, isotropic model by choosing a specific 
form of the quadratic free elastic energy density. 

The question of boundary conditions in higher gradient elasticity models has been a subject of constant 
attention. The matter is that in a higher gradient model, it is not possible to independently vary the test 
function and its gradient. Some sort of split into tangential and normal parts is usually performed (see e.g. 
EH EH EH). This is well known in general higher gradient models. The boundary conditions in the general 
case of gradient elasticity and strain gradient elasticity have been settled in the paper by Bleustein [5] , see also 
mrnrnm- However, as it turns out, the boundary conditions obtained by Tiersten and Bleustein in [85.! 
with respect to the special case of the indeterminate couple stress model are not the only possible ones in the 
framework of the indeterminate couple-stress model. 

While the strain gradient framework necessitates to work with a third order hyperstress tensor, the inde¬ 
terminate couple stress model is apparently simpler: it restricts the form of the curvature energy and allows 
to work with a second order couple-stress tensor work-conjugate to gradients of rotation. For this apparent 
simplification the indeterminate couple stress model has been heavily investigated and is still being heavily 
used as well. A first answer as regards boundary conditions has been given by Mindlin and Tiersten as well 
as Koiter EH 37] who established (correctly) that only 5 geometric and 5 traction boundary conditions can be 
prescribed. Their format of boundary conditions has become the commonly accepted one for the couple stress 
model EH EH S3 IH EH EH, all these papers using the same set of (incomplete) boundary conditions. It seems, 
to us, however, that the state of the art in general strain gradient theories [59l [86l SZJ E21 ESI [61 [561 EH SH EH 
is much more advanced as far as boundary conditions are concerned. 

This paper has been motivated by our reading of [32] 1U3 SI] EH- 25 321, in which the form of traction 
boundary conditions in the indeterminate couple stress model, together with an apparently plausible physical 
postulate lead to unacceptable conclusions, see m- Therefore, there had to be an underlying problem which 
we believe to have tracked down to the hitherto accepted format of boundary conditions. 

The main result of this paper, consisting in setting up a "stronlgy independent" set of boundary conditions 
for the couple stress model, has been announced in |66| . This contribution is now structured as follows: after a 
subsection fixing the notations used throughout the paper, we outline some related models in isotropic second 
gradient elasticity and we give a brief digression concerning differential geometry. In Section [2] we present the 
equilibrium equations and the constitutive equations of the indeterminate couple stress model as they have been 
derived in the literature. We also present the classical "weakly independent" boundary conditions proposed by 
Mindlin and Tiersten [59] and the main arguments of their proposal. Since we remark that these boundary 
conditions are not the only possible ones, in Section [3] we obtain the novel set of boundary conditions in the 
indeterminate couple stress model. To this main aim of our paper, we follow two different paths. On the 
one hand, we consider the indeterminate couple stress model as a special case of the second gradient elasticity 
model and we derive the "strongly independent" boundary conditions which follow naturally as restriction of 
such general framework, see Subsection |3.2| This kind of approach involves the third order hyperstress tensor 
as a reminiscence of the second gradient elasticity approach. However, in Subsection |3.3[ we prove that the 
equilibrium equations and the boundary conditions may all be rewritten in terms of Mindlin’s second order 
couple stress tensor. On the other hand, following the line of Mindlin’s argument in combination with some 
calculations specific to second gradient elasticity model, we set up a "direct approach" which leads to the 
same set of boundary conditions with those coming from second gradient elasticity, see Section |3.4| However, 
these boundary conditions do not always coincide with those proposed by Mindlin and Tiersten and which are 
accepted and used until now in the literature. We explain this fact in Subsection |1,4| where we explicitly show 
that, if an "a priori" equivalence can be found in most cases between our approach and Mindlin’s one, this is 
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not the case when considering "mixed" boundary conditions, simultaneously assigning the force and the curl 
of displacement (Mindlin) or the force and the normal derivative of displacement (our approach) on the same 
portion of the boundary. 

In the appendix, we give some explicit or alternative calculations which are used in the main text and we an¬ 
swer again to the question: what are the missing steps in Mindlin and Tiersten’s approach? (see Appendix | A. 4[ ). 
We do so by using all the arguments provided throughout the paper and pointed out in different circumstances. 
We end our paper by some concluding diagrams summarizing our findings. 


1.1 Notational agreements 

In this paper, we denote by R 3x3 the set of real 3x3 second order tensors, which will be written with 
capital letters. We denote respectively by • , : and (•,•) a simple and double contraction and the scalar 
product between two tensors of any suitable ordeiQ Everywhere we adopt the Einstein convention of sum 
over repeated indices if not differently specified. The standard Euclidean scalar product on R 3x3 is given by 
(A, y) R 3x3 = tr(X -Y T ), and thus the Frobenius tensor norm is ||X|| 2 = (X, X) R 3x3. In the following we 
omit the index R 3 ,R 3x3 . The identity tensor on R 3x3 will be denoted by 1. so that tr(X) = (A", 1). We 
adopt the usual abbreviations of Lie-algebra theory, i.e., so(3) := {A £ R 3x3 |X T = —A} is the Lie-algebra 
of skew symmetric tensors and sl(3) := {A' £ R 3x3 | tr(A) = 0} is the Lie-algebra of traceless tensors. For 
all A' £ R 3x3 we set symX = \{X T + X) £ Sym, skew A = |(A — X T ) £ so(3) and the deviatoric part 
dev A = A — | tr(X) 1 £ sl(3) and we have the orthogonal Cartan-decomposition of the Lie-algebra fll(3) 

g[(3) = {sl(3) fl Sym(3)} ® so(3) ® R-l, A = devsymX + skew A + ^ tr(X) 1. (1.1) 

Throughout this paper (when we do not specify else) Latin subscripts take the values 1, 2,3. Typical conventions 
for differential operations are implied such as comma followed by a subscript to denote the partial derivative 
with respect to the corresponding cartesian coordinate. Here, for 

( 0 -a 3 a 2 \ 

a 3 0 —ai I £ so(3) (1-2) 

-a 2 ai 0 / 

we consider the operators axl : so(3) —> R 3 and anti : R 3 —> so(3) which verify the following identities 

axl(A) := (ai, a 2 , a 3 ) T , A. v = (axl A) x v, (anti =—eij k v k , Vv£R 3 , (1-3) 

(axl A) k = -^ CijkAij = i tkijAji , Aij = -e ijk (axl A) k =: anti(axl A)^, (a x b)i = e ijk ajb k , 

where e ljk is the totally antisymmetric Levi-Civita third order permutation tensor. 

We consider a body which occupies a bounded open set fl of the three-dimensional Euclidian space R 3 and 
assume that its boundary <9f1 is a smooth surface of class C 2 . An elastic material fills the domain fl C R 3 and 
we refer the motion of the body to rectangular axes Oxi. 

With reference to FigjTJ n is the outward unit normal to <9fl, T is an open subset of the boundary di 1, v~ 
is a vector tangential to the surface dil \ T and which is orthogonal to its boundary d(dfl \ T), t~ = n x v~ is 
the tangent to the curve d(dLl \ T) with respect to the orientation on dfl \ T given by the outward unit normal 
n to this surface. Similarly, v + is a vector tangential to the surface T and which is orthogonal to its boundary 
9T, t + = n x v + is the tangent to the curve dT with respect to the orientation on T. 

In the following, given any vector field a defined on the boundary <9f2 we will also set 

[(«, v )\ := {a + ,v + ) + (a~,v~) = (a + ,v) - {a~,u) = (a + - a~,v), (1.4) 


which defines a measure of the jump of a through the line 9T, where v := v + = — v and 


lim 

ie3!i\r 
x —y dr 


:= lim [•]. 
x e r 
or 


1 For example, (A ■ v)i = AijVj , (A • B )^ 
v • w = ViWi, (A, B) = AijBij etc. 


AijBjk, A B — AijBji, (C • B^ijk — C ijp B ph , (C : B^i — C^jpBpj , — 
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Figure 1: The domain 12 C M 3 together with the part T C dfl, where Dirichlet boundary 
conditions are prescribed. We need to represent the boundary conditions on a disjoint union 
of dfl = (dfl \ T) U T U ar, where T is a open subset of dfl. 


We are assuming here that dfl is a smooth surface. Hence, there are no geometric singularities of the boundary. 
The jump [•] arises only as a consequence of possible discontinuities which follows from the prescribed boundary 
conditions on T and dfl \ T. Nevertheless, if one would like to explicitly consider continua with non-smooth 
boundaries, the jump conditions to be imposed at the edges of the boundary would be formally the same to 
those that we will present in the remainder of this paper, with the precision that the jump would in this case 
indicate a true jump across a geometrical discontinuity of the surface. 

The usual Lebesgue spaces of square integrable functions, vector or tensor fields on fl with values in R, R 3 
or R 3x3 , respectively will be denoted by L 2 (fl). Moreover, we introduce the standard Sobolev spaces [21321 IS] 


H 1 (f2) = {u £ L 2 (fl) | grad it € . 
H(curl; H) = {v S L 2 (fl) | curlu 
H(div; fl) = {v £ L 2 (fl) | div v £ 


^)}, 

M Hqn) := 

£L 2 m, 

M H(curl;fi) 

L 2 m, 

M H(div ; n) : 


u lll 2 (Q) + lls rac Ml! 2 (0) > 

= IMIi 2 (fi) + IIcurl^||, (1-5) 

: IMIi 2 (Q) + 11 div11 ^2 (f2 ), 


of functions u or vector fields v, respectively. 

For vector fields v with components in i.e. v = (iq, u 2 , u 3 ) T ,Vi £ H^fl), we define 

Vt = ((V V\) T , (V t> 2 ) T , (V ii 3 ) t ) , while for tensor fields P with rows in H(curl;fl), resp. H(div;fi), i.e. 

P = (Pi , P 2 , Pj)'. Pi £ H(curl; fl) resp. Pi £ H(div ; fl) we define CurlP = ((curl Pi) T , (curlP 2 ) T , (curlP 3 ) T ) T , 
DivP = (div P ll div P 2 , divP 3 ) T . The corresponding Sobolev-spaces will be denoted by 


H 1 (S7), H 3 (Div ; fl), H^Curl; H). 


1.2 The indeterminate couple stress model 

In the indeterminate couple stress model we consider that the elastic energy is given in the form 

W = Hji n (Vu) + W C urv(V[axl(skew Vu)]) = Hji n (Vu) + IT curv (Vcurlu), (1.6) 


where 


\ K 

Hji n (Vu) = [i || sym Vu|| 2 + — [tr(sym Vu)] 2 = p || devsym Vtt|| 2 + — [tr(sym Vu)] 2 . (1.7) 

Here, /i > 0 is the infinitesimal shear modulus, k = 2 ^i( 3A > 0 is the infinitesimal bulk modulus with A the 
first Lame constant. In order to discuss the form of the curvature energy W curv (V[axl(skew Vu)]), let us recall 
some variants of the linear isotropic indeterminate couple stress models. Some parts of this classification have 
already been included in the paper m but we include them also here for the sake of completeness: 

• the indeterminate couple stress model (Grioli-Koiter-Mindlin-Toupin model) |55I [21 371 [551 [571 15^1136) 
in which the higher derivatives (apparently) appear only through derivatives of the infinitesimal continuum 
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rotation curl it. Hence, the curvature energy has the equivalent forms 


G!i 

Wcurv(V[axl(skew Viz)]) = — || symVcurl it|| 


012 


skewVcurl it || 2 


= ol\ || sym V[axl(skew Vu )]|| 2 + a 2 || skew V[axl(skew Vit)]||" 


OL\ < '~*-z 

= — || dev symVcurl tx || 2 + || skewVcurl u 


0L2 


( 1 . 8 ) 


Here, we have used the identities 
2 axl(skewVit) = curl it, 


tr[V[axl(skew Vu)]| = - tr[V[curlu]] = -div[curlu] = 0, (1.9) 


together with the fact that Vcurl it is a trace-free second order tensor and hence so is symVcurl u. This 


implies that devsym Vcurlu = symVcurl it. Although this energy admits the equivalent forms (1.8 11 and 


( 1 . 813 , the equations and the boundary value problem of the indeterminate couple stress model is usually 


formulated only using the form (I. 8 I 1 of the energy. Hence, we may individuate one of the aims of the 


present paper in the fact that we want to formulate the boundary value problem for the indeterminate 
couple stress model using the alternative form (1.8 12 of the energy of the Grioli-Koiter-Mindlin- 


Toupin model (see Section [ 2 ]) . We also remark that the spherical part of the couple stress tensor is zero 
since tr(Vcurlu) = div(curlit) = 0. In order to prove the pointwise uniform positive definiteness it is 
assumed, following 03 . that a 1 > 0 ,02 > 0 (corresponds to — 1 < 77 := < 1 in the notation of |47|). 

Note that pointwise uniform positivity is often assumed m when deriving analytical solutions for simple 
boundary value problems because it allows to invert the couple stress-curvature relation. We have shown 
elsewhere [3D] that pointwise positive definiteness is not necessary for well-posedness. 

• In this setting, Grioli [351 1.16] (see also Fleck |2511261I27| ) initially considered only the choice oq = 0 . 2 - In 
fact, the energy originally proposed by Grioli m is 

r 2 “1 


W CU rv(D 2 u) = nL 2 c — [ ||V(curl it )|| 2 + r/ tr[(V(curl u)) 2 }] 

= nL 2 ^p[|| devsym V[axl(skewVu )]|| 2 + || skew V[axl(skew Vit)]|| 


= fiL 


2 Q, i 
4 L 


- 77 (V[axl(skew Vit)], (V[axl(skew Vit)]) T )] 
dev symV(curl u )|| 2 + || skew V(curl u )|| 2 

+ 77 (V(curl u), (V(curl u)) T ) 


( 1 . 10 ) 


= fiL 2 -^- [(1 + 17 ) || devsym V(curl u )|| 2 + (1 — 77 ) || skew V(curl u) 


Mindlin [5D] p. 425] (with 77 = 0) explained the relations between Toupin’s constitutive equations [ 5P] 
and Grioli’s [ 35] constitutive equations and concluded that the obtained equations in the linearized theory 
are identical, since the extra constitutive parameter 77 of Grioli’s model does not explicitly appear in the 
equations of motion but enters only the boundary conditions, since V axl(skewVu) = [Curl (sym Vit)] T , 
Div Curl (•)= 0, and 

Div{antiDiv[Vaxl(skew Vzt)] T } = Div{antiDiv[Curl (sym Vu)]} = Div{anti(0)} = 0. 


The same extra constitutive coefficient appears in Mindlin and Eshel m and Grioli’s version (|1.10[). 


• the modified - symmetric couple stress model - the conformal model. On the other hand, in the 
conformal case ESI CD one may consider that 02 = 0 , which makes the second order couple stress tensor 
m symmetric and trace free m- This conformal curvature case has been considered by Neff in m, the 
curvature energy having the form 


Wcurv(Vcurl It) 


Ol 

4 


|| symVcurl it|| 2 . 


( 1 . 11 ) 
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Indeed, there are two major reasons uncovered in m for using the modified couple stress model. First, in 
order to avoid singular stiffening behaviour for smaller and smaller samples in bending [Bjl] one has to take 
a 2 = 0. Second, based on a homogenization procedure invoking an intuitively appealing natural “micro¬ 
randomness" assumption (a strong statement of microstructural isotropy) requires conformal invariance, 
which is again equivalent to 02 = 0. Such a model is still well-posed 051 leading to existence and 
uniqueness results with only one additional material length scale parameter, while it is not pointwise 
uniformly positive definite. The initial motivation of Yang et al. [ 88 ] for using the modified couple stress 
model is based on incomplete arguments m : even if their conclusions concerning a symmetric couple 
stress tensor may be kept in some particular phenomenological cases m- 

1.3 Brief digression concerning differential geometry 

When dealing with higher order theories it is suitable to introduce (see also [81 j ;T5] for details) two second order 
tensors T and Q which are the two projectors on the tangent plane and on the normal to the considered surface, 
respectively. As it is well known from differential geometry, such projectors actually allow to split a given vector 
or tensor field in one part projected on the plane tangent to the considered surface and one projected on the 
normal to such surface. We can introduce the quoted projectors as 

T = t®t + v®v = \—n®n, Q = n®n. 

It is easy to check that the following identities are verified by the two introduced projectors 

T + Q = 1 , T -T = T, QQ = Q, T-Q = 0, T = T t , Q = Q T . (1.12) 

It is then straightforward that any vector v can be represented in the local basis {r, n, ix} as 

v = v ■ (T + Q) = (T + Q) ■ v = (v, t) t + (v , v) v + (v, n) n , 

T-v Qv 


or equivalently the components of v can be written as 


Vi — Vh Thi Vh Qhi — Vh T Qih Vh- 

Analogously, a second order tensor B can be represented in such local coordinates a s B = (T + Q) T ■ B ■ (T + Q) 
and this representation can be also straightforwardly generalized to a tensor of generic order N. To the sake 
of generality, from now on we introduce a global orthonormal basis {e 1 , e 2 , e 3 } in which all fields (included r, v 
and n) will be represented if not differently specified. This also implies that all the space differential operations 
which will be performed in the following calculations are all referred to the coordinates (Xi, X 2 , X 3 ) associated 
to such global basis. 

From m p. 58, ex. 7], we have the following variant of the surface divergence theorem 
Proposition 1.1. Suppose that w satisfies (u>,n) = 0 on the surface S C R 3 , then 

/ (T,'\7w)da= / (■ nxw,r)ds , (1.13) 

Js JdS 

where t is a vector tangent to S and to dS. 

Taking now w = T ■ v, for arbitrary v £ R 3 and using that 

(n x (T ■ v ), r) = (r x n,T ■ v) =: (u, T ■ v) = (T ■ v, v) 

we obtain 

[ (T, V (T • v)) da = f (n x (T ■ v),t) ds = f (T ■ v, v) ds = f (v,v)ds, (1.14) 

JS JdS JdS JdS 

for v £ R 3 , since T ■ v = v ■ T = v. We explicitly remark that the vector v is tangent to the surface S but 
orthogonal to its boundary dS (it points inward or outwards, depending on the choice of the orientation of n 
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and t: usually, n is chosen to be the outward unit normal to S and the orientation of r is chosen in order to 
have that the vector u, orthogonal to the border of S, is pointing outward the border of the surface S itself). 
With such notations, when considering a vector field v defined in the vicinity of the considered surface, the 
surface divergence theorem can be applied to the projection of v on the tangent plane to the surface as follows 
(see e.g. [55]) 

f Div s (v) da := f (T, V (T ■ v)) da = f (T-v,v)da= f {v,T ■ v) da = f (v,v)ds , (1.15) 

Js Js JdS JdS JdS 

where, clearly, in the last equality the notation 

Div s (v ) := (V (T-v),T) 

for the surface divergence of a generic vector v has been used. Equivalently, in index notation the surface 
divergence theorem reads 

/ (TijVj) k T ik da = / ViT ik v k ds= / ViVids. (1.16) 

JS ’ JdS JdS 

This definition of surface divergence can be extended to higher order tensors, in particular, for a second order 
tensor B , its surface divergence is introduced as the vector of components [Div s (13)] i = (T • B)ij k Tj k . 

Remark 1.2. We explicitly remark that if S coincides with the boundary dfl of the considered body and T is 
an open subset ofdfl, then the surface divergence theorem (1.151 implies (see Fig. [7] and eq. ( |1.4 | ) ) 

I Di v s (v) da := f (T, V (T • v)) da = f Div s (u _ ) da+ I Div s (u + ) da (1.17) 

Jd n Jd n JdQ\ r Jr 


') ds 


Jd(d n\r) 

Clearly, we can equivalently write in index notation 


/an\r 

/ {V' , v 
Id r 


+ ,,+ '' ds = 


/ 
Jd r 


I ds. 


/ (TijVj) k Tik da = [ Vj Vj\ ds. 

Jd n ’ Jd r 


(1.18) 


1.4 Weakly and strongly independent surface fields 

Since it will be useful in the following, we give in this subsection some definitions which will be used throughout 
the paper. In particular, we introduce the notions of "strongly" and "weakly independent" vector fields defined 
on suitable regular surfaces. 

Definition 1.3. [Weakly independent surface fields] Given two vector fields u and v defined on a suitable regular 
surface S, we say that they are "weakly independent" if we can arbitrarily assign u and v, independent of each 
other, by choosing two vectors u and v and a function f such that 

u = u and v = f(u,v). (1.19) 

By means of the notion of weak independence, we can arbitrarily fix the vectors u and v on the surface S, but a 
variation of the chosen u induces a variation on one part of v. Nevertheless, the vector v can still be arbitrarily 
chosen thanks to the arbitrariness ofv and f. 

Definition 1.4. [Strongly independent surface fields] Given two vector fields u and v defined on a suitable 
regular surface S, we say that they are "strongly independent" if we can arbitrarily assign u and v, independent 
of each other, by choosing two vectors u and v such that 

u = u and v = v. (1.20) 

The notion of strong independence allows to arbitrarily fix the two vector fields u and v on the surface S and a 
variation in the choice of the first vector u = u does not induce a variation in the vector v. 



Considering the power of external actions for a second gradient continuum 

V ext = f ( t ext ,6u)+ f ( g ext ,D\Su )), 

Jon Jon 

where Dqciu) is a suitable first order space differential operator, we say, by extension of the former definition, 
that t ext and g ext are strongly and weakly independent tractions whenever they are the conjugates of strongly 
or weakly independent surface vector fields Su and D 1 (Su). 

An example of strong and weak independence which will be pertinent in the framework of the present paper 
is that which can be established between the displacement assigned on the surface and its curl or, alternatively, 
its normal derivative. Indeed, if for convenience we choose an orthonormal local basis {r, u, n} on S, where r 
and v are unit tangent vectors to S, while n is the outward unit normal vector, we can recognize that 


Vw • n = 


(du T du v du r 


\dx n ’ dx n ' dx r 


and curl u = 


/ du n 

du v 

du T 

du n 

du„ 

du T \ 

ydx,. 

dXn 

dx n 

dx T ’ 

dx T 

dx v ) 


( 1 . 21 ) 


where we clearly indicated by u T , u u and u n the components of the displacement field in such local basis and 
x T , x v and x n the space coordinates of a generic material point in the same reference frame. 

It is known that the fact of fixing the displacement field u on the surface does not fix its normal derivatives 
du T /dx n , du v /dx n , du n /dx n , while it fixes the tangential derivatives du n /dx u , du n /dx T1 du„/dx T , du T /dx v . 
In this optic, we can say that u and V« • n are strongly independent surface vector fields, while u and curlu 
are weakly independent surface vector fields. Indeed, even if when fixing the displacement u on the surface its 
tangential derivatives result to be fixed, the arbitrariness of the normal derivatives still allows to globally choose 
curl u in an arbitrary way. 

Therefore (regarding the curlu), it is impossible to fix as constant the curlu but varying only u. In contrast, 
it is possible to fix as constant S7u.n and to vary u arbitrarily. In summary in this example 


u and Vtt. n 


are strongly independent, 


u and curl u 


are (only) weakly independent. 


2 Bulk equations and (Mindlin’s weakly independent) boundary con¬ 
ditions in the indeterminate couple stress model 

In this section we re-propose the results presented by Mindlin and Tiersten but explicitly using a variational 
procedure. More particularly, we try to explicitly present some reasonings that can be made to obtain their 
bulk equations and boundary conditions. Nevertheless, how we will show in due detail in the remainder of this 
paper, such reasonings are not in complete agreement with other sets of boundary conditions which can be 
assigned when dealing with a couple stress theory, and which are equally legitimate. 

2.1 Equilibrium and constitutive equations of the indeterminate couple stress 
model 

Since we consider that the solution belongs to H 1 (il), we take free variations Su £ of the energy 

W{'\7u 1 Vcurlu) = Mdi n (Vrt) + Wc Urv (Vcurlu), (2-1) 


where 


Wi in (\7u) =n || symVu|| 2 + ^ [tr(Vu)] 2 = g || devsymVu|| 2 + 2 ^ ^ 3 - [tr(Vu)] 2 , 
Wcurv(Vcurlu) =a\ || devsym V[axl(skew Vu)]|| 2 + «2 || skew V[axl(skew Vtt)]|| 2 , (2.2) 
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and we obtain the first variation of the action functional as 
SA = S / — W(Vu, Vaxl(skewVu)) dv = — / 2 / u(symVu,symVi5u) + Atr(Vu) tr(V<5u) 

J Q J i 1 . 

+ 2ai [(devsymV[axl(skew Vu)],devsymV[axl(skew Vdu)]) (2.3) 
+ 2a 2 (skew V[axl(skew Vu)],skew V[axl(skew V<5u)])] dv. 

Or equivalently, applying the classical divergence theorem 

SA = / (Div (2 p, sym Vu + Atr(Vu)l), 8u) dv — / ((2/r sym V« + Atr(Vu)l) • n, Su) da 
J n Jan 


-2 


oq devsym[V axl(skew Vu)] + a 2 skew[V axl(skew Vu)], V axl(skew V<5u)]y dv. (2.4) 


The classical divergence theorem can be applied again to the last term appearing in the previous equation. In 
particular, we can notice that the following chain of equalities holds 


aq devsym[V axl(skew Vu)] + a 2 skew[V axl(skew Vu)], V axl(skew Vdu)]y dv 

— Div ( oq dev sym[V axl(skew Vu)] + a 2 skew[V axl(skew Vu)]j , axl(skew V<5u)] ^ dv 


1 an 


aq devsym[V axl(skew Vu)] + a 2 skew[V axl(skew Vu)] 


n, axl(skew V<5u)y da 


— anti{Div(oq devsym[V axl(skew Vu)] + a 2 skew[V axl(skew Vu)])}, skew Vduy dv 


1 an 


aq devsym[V axl(skew Vu)] + a 2 skew[V axl(skew Vu)] 


n, axl[skew V<5u] ) da 


— skewanti{Div(«i devsym[V axl(skew Vu)] + a 2 skew[V axl(skew Vu)])}, Vduy dv 


+ 


lan 


aq devsym[V axl(skew Vu)] + a 2 skew[V axl(skew Vu)] 


n, axl(skew V<5u)y da 


Div [skewanti{Div(ai devsym[V axl(skew Vu)] + a 2 skew[V axl(skew Vu)])}], SuJ dv 
j [skew anti{Div(ai dev sym[V axl(skew Vu)] + a 2 skew[V axl(skew Vu)])}] • n, Su ^ da 


lan 


ct\ devsym[V axl(skew Vu)] + a 2 skew[V axl(skew Vu)] 


• n, axl(skew Vdu)y da, 


where n is the unit outward normal vector at the surface dil. Hence, recalling that skewanti(-) = anti(-), the 
first variation (2.4) of the action functional can be finally rewritten as 

SA= / (Div (er — r), Su) dv — / ((a — r) • n,5u) dv — / (m • n, axl(skew Vdu)) da, 

Jn Jan Jan 


for all virtual displacements Su £ C°°(f2), where a is the symmetric local force-stress tensor 

er = 2 /r sym Vu + A tr(Vu)l £ Sym(3), 


(2.5) 
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r represents the second order nonlocal force-stress tensor (which here is automatically skew-symmetric) 

t = a i antiDiv(devsym[V axl(skew Vu)]) + a 2 antiDiv(skew[V axl(skew Vu)] (2-6) 

= — antiDiv(devsym[Vcurly]) + — antiDiv(skew[Vcurly]) 


GL 2 

— devsym(Vcurly) + — skew(Vcurly) 


= anti Div 
= ^ anti Div [m] £ so(3), 


and 


m = 2 ot\ dev sym( V axl(skew Vu)) + 2 02 skew(V axl(skew Vu)) 

= [ai devsym(Vcurlu) + a 2 skew(Vcurlu)] (2.7) 

<*1 + 0:2 _ , ai — a 2 , ,t 

= --- Vcurl u H--- [Vcuri u\ 

is the second order hyperstress tensor (couple stress) which may or may not be symmetric, depending on the 
material parameters. The asymmetry of force stress is a hidden constitutive assumption, compared to the 
development in [30 ]. Postulating a particular form of the power of external actions, the equilibrium equation 
can therefore be written as 


Div ototai + / ext = 0, 


( 2 . 8 ) 


where we clearly set <7 t otai = cr — r. 


2.2 Classical (weakly independent) boundary conditions in the indeterminate cou¬ 
ple stress model 

In the previous section, we have shown that the power of internal actions can be finally written in compact form 
as 

r p int = SA= f (Div(a — r), Su) dv — f ((a — r) ■ n, Su) da — ^ f (rh ■ n, curl Su) da. (2.9) 

J n Jon 2 J dn 

Such form of the power of internal actions seems to suggest 6 possible independent prescriptions of mechanical 
boundary conditions; three for the normal components of the total force stress (<7 — t) • n and three for the 
normal components of the couple stress tensor. The possible Dirichlet boundary conditions on T C dfl would 
seem to be the 6 condition^ 

u = zt ext , axl(skewVu) = -curl u = a ext (or equivalently curlw = 2 a ext ), ( 2 . 10 ) 


for two given functions it ext ,a ext : R 2 3 —>• R 3 at the boundary T C dfl (3+3 boundary conditions). 
However, following Koiter m we must note that the following remark holds true: 


Remark 2 . 1 . Assume that u £ C , °°(f2) and u 


is known. Then curly 


exists and for all open subsets T C dfl 


the integral J r {curlw, n) da is already known, while f r (curlu,r) da is still free, where r is any tangential vector 
field on T C dVL. This fact follows using Stokes’ circulation theorem 


/ (curly, n) da = / (u('y(t)),'y'(t))ds, 
Jr Jd r 


( 2 . 11 ) 


where t is a continuous unit vector field tangent to the curve <9T = {7 (t) |f £ [a, b]} compatible with the unit 
vector field n normal to the surface T. 


This leads us to the next correct observation: 

2 as indeed proposed by Grioli m in concordance with the Cosserat kinematics for independent fields of displacements and 

microrotation. 
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Remark 2.2. Only the two tangential components o/curlrt may be independently prescribed on an open subset 
of the boundary. However, we may have six independent conditions in one point on T, but not on an open subset 
of the boundary. 


Already Mindlin and Tiersten |59| have also correctly remarked that in this formulation only 5 mechanical 
boundary conditions can be prescribed. Using our notations, their argument runs as follows: 

-(to • n, curlfe) = — ((T + Q) ■ (fh ■ n),curl<5«) = -(T- (to - n), curl for) + -(Q • (to • n), curlew) 

= —{T ■ (to ■ n), curlJu) + - ((n <S> n) ■ ( fh ■ n ), curl Su) = —{T ■ (fh ■ n), curlew) + -((((symm) • n, n)) n, curlckt) 

= — (T ■ (to • n), curlJu) + -(n, (((symTO.) • n, n)) curl Su) (2.12) 

= —{T ■ (fh • n), curl Su) + — (n, curl [(((symm) • n, n)) Su]} — -(n, V (((symm) • n, n )) x 5u) 

= -(T • (to • n), curlJu) + — (n, curl [(((symm) • n, n)) <5«]) — -(n x V (((symm) • n, n )), <5it) 

= -(T ■ (fh ■ n), curlew) + — (n, curl [(((symm) • n, n)) <5u]) — -(n x V (((symm) • n, n )), Su), 

where ® denotes the dyadic product of two vectors, we have used the property (77 <g> £) • a = rj (£, a) (for vectors 
77 , ^ and a), the formula curl (ip Su) = \7ip x Su + ip curl i5u (for any scalar field ip ) and the fact that n®n is a 
symmetric second order tensor. The power of internal actions (2.91 can hence be rewritten as 


pmt _ / ^j3iv(,T — 7 -), Su) dv — / ((a — t) ■ n, Su) da 


/an 


f | -(T ■ (to • n), cxvrlSu) + -(n, curl [(((symm) • n,n)) (5u]) — — (n x V (((symm) • n,n)) ,Su) 1 da. 
Jon l 2 2 2 J 


It can be remarked that the second term in the last surface integral can be rewritten as a bulk integral by means 
of the divergence theorem, so that the power of internal actions can also be rewritten in a further equivalent 


form in addition to the one already established in (2.9) 


V int = / (Di v(a — t),5u) dv — / - div {curl [(((symm) • n, n)) <5w]} dv 
J n J n 2 

— f ((a — t) ■ n — -n x V (((symm) • n, n)) , Su) da — j -((m • n), T ■ curl Su) 
Jd n 2 Jd n 2 


(2.13) 


da, 


where the fact that the tangential projector T is symmetric has also be used. 

Mindlin and Tiersten |59| concluded that 3 boundary conditions derive from the first surface integral and 
two other from the second surface integral, since [59[ p. 432] “the normal component of the couple stress vector 
[(to • n,n) = ((symfh.) • n,n)\ on dfl enters only in the combination with the force-stress vector shown in the 
coefficient of Su in the surface integral ...”. Indeed, Mindlin and Tiersten are assuming to assign arbitrarily the 
displacement and the tangential components of its curl on the surface dfl. 

As we will deeply discuss in Section [4] this choice leads to a possible set of boundary conditions in the 
indeterminate couple stress model. Nevertheless, this choice is not unique and assigning at the boundary 
different virtual fields as the virtual displacement and its normal derivative will lead us to a set of boundary 
conditions that are not equivalent to those proposed by Mindlin and Tiersten. 


2.2.1 Geometric (essential, or kinematical) weakly independent boundary conditions 

Based on the expression ( |2.13 1 of the power of internal actions, Mindlin and Tiersten [59] concluded that the 
geometric boundary conditions on T C dfl are the five independent conditions 


= u 


(1 — n ®n) ■ curlitL = (1 — n®n)-a e 


(3bc) 

(2bc) 


(2.14) 
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for given functions u ext , a ext : R 3 —>• R 3 on the portion T of the boundary. 

An equivalent form of the above boundary condition is 

uL = u ext , (3bc) 

i t (2.15) 

(1 — n<8> n) ■ axl(skew Vu)| r = |(1 — n® n) ■ a ext , ( 2 bc) 

for given functions u ext ,a ext : R 3 —> R 3 at the boundary. The latter condition prescribes only the tangential 
component of axl(skewVu) = -( curly. Therefore, we may prescribe only 3+2 independent geometric boundary 
conditions. Regarding this formulation, an existence result was proven in [45] . 

In order to give a first comparison with the boundary conditions which are coming from the full strain 
gradient approach, let us remark that: 

Lemma 2.3. [Equivalence of geometric boundary conditions] We consider a vector field 
u : R 3 —> R 3 , u £ C°°(£l). The following sets of boundary conditions are equivalent: 


(1 — n < 8 > n) ■ curl u 

in the sense that one set of boundary conditions defines completely the other set of boundary conditions, where 
n is the unit outward normal vector at the surface T C dfl and a ext and a ext can be a priori related. 

Proof. The proof is included in Appendix |A.5| □ 

In summary, we can say that if purely kinematical boundary conditions are assigned in the indeterminate 
couple stress model, in virtue of the previous Lemma |2.3[ it is equivalent to assign on one portion of the 
boundary the displacement and the tangential part of its curl or the displacement and the tangential part of its 
normal derivative. As we will see, things become much more complicated when one wants to deal with traction 
or mixed boundary conditions, since it is not straightforward to individuate the equivalence between different 
sets of boundary conditions which are nevertheless equally legitimate. 


= ( 1 - 


+> 


n <+ n) - a 


(1 — n (g> n) • Vu ■ 


= ( 1 - 


n <+ n -a 


(2.16) 


2.2.2 Classical (weakly independent) traction boundary conditions 


Always considering the expression (2.131 of the power of internal actions, the possible traction boundary con¬ 
ditions on dfl \ r given by Mindlin and Tiersten m in the axial formulation are 


(a — t) ■ n— -n x V (((sym m) ■ n,n )) = f ext , 

(1 — n ® n) • rh ■ n = (1 — n 0 n) ■ g ext , 


traction 

double force traction 


(3 be) 

(2 be) (2.17) 


for prescribed functions t ext , g ext : R 3 —> R 3 on the portion dfl \ T of the boundary. 

Since Su and (1 — n ® n) ■ curl Su are weakly independent, at this point we are tempted to conclude that the 
equality 


/ ( t ., Su) da + / (g, (1 — n (g) n) ■ curl Su) da = 0 

ldQ\T ' JdQ\T 


(2.18) 


for all Su £ C 2 (n) does not imply that t 
using the Lemmas included in Appendix 


dfi\r 

O] 


= 0 and (1 — n®n) g 


an\v 


= 0. However, this holds true after 


We want also to explicitely point out that (2.141 and (2.17) correctly describe the maximal number of 


independent boundary conditions in the indeterminate couple stress model but even if these conditions have 
been re-derived again and again by Yang et al. [ 88 ], Park and Gao m , [50| . etc. among others they are not the 
only possible choice in the couple stress model. This is explained in the following two subsections. Prescribing 
Su and (1 — n ® n) ■ curl Su on the boundary means that we have prescribed independent geometrical boundary 
conditions, this is also the argumentation of Mindlin and Tiersten [52[, Koiter m, Sokolowski[82j, etc. However, 
the prescribed traction conditions are not stronlgy independent but only weakly independent in the sense 
established in Section [+4] For this reason we claim that, in order to prescribe strongly independent geometric 
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boundary conditions and their corresponding energetic conjugate, we have to prescribe u and (l-n®n)-Vu-n. 
In other words, as it is well assessed in the framework of full second gradient theories (see also the expression 
of the external power given in (3.8)), we prescribe on the boundary the following part of the power of external 
actions 


i an 


(f ext , Su) da + / {g cxt , (1 — n <8> n) ■ V5u ■ n) da, 


i an 


(2.19) 


in which now Su and (1 — n ® n) ■ Vdu ■ n are strongly independent and <? ext does not produce (anymore hidden) 
work against Su. This type of strongly independent boundary conditions are also correctly considered already 
by Bleustein j9], but for the full strain gradient elasticity case only. We give more details in the following section. 


3 Through second gradient elasticity towards the indeterminate cou¬ 
ple stress theory: a direct approach 

Independently of the method that one wants to choose to set up the correct set of bulk equations and associated 
boundary conditions for the indeterminate couple stress model, such set of equations must be compatible with 


a variational principle based on the form (2.11 of the strain energy density. We present two different ways of 


performing such variational treatment: the first one passes through a full second gradient approach and the 
second one, which we call direct approach, is based on the fact that the curvature energy is regarded as a 
function of the second order tensor Vcurlu instead than of the third order tensor VVzt or V syrri V'u. 

Worthless to say, as expected, we will find that such two approaches are equivalent and we will explicitly 


establish this equivalence in eqs. (3.31)-(3.341. 


3.1 Second gradient model: general variational setting 

In this section, we show how the couple stress model can be regarded as a particular case of the second gradient 
model. 


3.1.1 First variation of the action functional: power of internal actions 

Let us consider the second gradient strain energy density W(Vu, VVu) and the associated action functional in 
the static case (no inertia considered here) 


A = - f W{SJu,S7S7u)dv. 
Jn 


The first variation of the action functional can be interpreted as the power of internal actions V lnt of the 
considered system and can be computed as follows 


V int = 5A=- 


/ dW 

\ dui^j 


Sui. 


dW 


du. 


i,jk 


Suijk ) dv, 


where we used Levi-Civita index notation together with Einstein notation of sum over repeated indices. Inte¬ 
grating a first time by parts and using the divergence theorem we get 


8A = - 


dW 


an ou i,j 


rijSui da 


dW 


du. 


Sui dv — 




r ow 

Ian 


rikSui j da + 


dW 


du. 


Sui j dv. 


Integrating again by parts the last bulk term we get 


8A = - 


i an 


dW 
du 


hJ 


dW 


du. 


bifc / ,fc 


rijSui da - 


dW 
du. 


i,3 


dW 


du. 


ijk J ,k 


Sut dv — 


i,jk / k 


f dW 

Ian duijk 


rikSui. 7 da, 


which can also be rewritten as 

SA = - 


/ (.<Tij - m ijk,k) rijSui da + (cr.y - rt \ijk,k) ,• Sui dv - / m ijk n k Suij da, (3.1) 

Jan Jn Jan 
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if one sets 


dW 
du 


rXijk — 


dW 


or equivalently, in compact notation: 




dW 


du 


i,jk 


dW 


m = 


dX/u SVVu' 

3.1.2 Surface integration by parts and independent variations 


(3.2) 


At this point, it must be considered that expression (3.1I can still be manipulated remarking that the tangential 
trace of the gradient of virtual displacement can be integrated by parts once again and that the surface divergence 
theorem can be applied to this tangential part of XSu. Using the brief digression concerning differential geometry 
and recalling the properties ( 1 . 12 ) of the tangential and normal projectors, we can now ulteriorly manipulate 
the last term in eq. (3.1) as follows, 

I ri k S'u>jj dct —. I BijSui^Shj du — I B^jSu^ ^ (I'h.y A Qhj) du 

J dQ J dQ J dQ 


I BijSui ^Thj da + I BijSui^Qhj da , 

J dQ J dQ, 

/ (T hj B ij )Su iyh da+ (BijUj) 6 (u ith ) n h du 

JdQ JdQ 

I (d^hpd^pjBij ) dui^h da + I {BijUj) S (ui^hUfa) da , 
JdQ JdQ 


(3.3) 


where we clearly set 


Bij — nr ijk rifc. 


(3.4) 


We can hence recognize in the last term of this formula that the virtual variation of the normal derivative 
Ui'h.nh =: it" = (Vtt • n)i of the displacement field appears. As for the first term, it can be still manipulated 
suitably integrating by parts and then using the surface divergence theorem ( 1 . 181 , so that we can write 


/an 


mp k. ///,: du, ! du 


/an 


’ hp 


ir^pjh (TpjBij) h Sui da + j (BijUj) Siii da , 


LpjJ^zjJ h 


/dQ 


(3.5) 


/ar 


lBi V v p 8ui\ds - / (T pj B i: j) h T hp Sui da + / (B^rij) Su™ da. 

JdQ ’ JdQ 


The final variation of the second gradient action functional given in (3.1 1 , can therefore be written as 
V int =SA= (<Tij - m ijk ,k) Sui dv - 


/an L 


/an 


(uij n Mjk,k) Tij (Tpjmij k n k ) ^ T kp 

(mij k n k n.j) Su^ da - / lxn i j k n k i'j8uijds, 

J ar 


Sil da 


(3.6) 


or equivalently in compact notation 

V mt = 5A = I (Div (a — Div m), Su) dv — f ((a — Div m) • n — V [(m • n) ■ T\ :T,8u) da 
in Jon 

— / ((m • n) ■ n, (X78u) ■ n) da — / [((m • n) ■ v, <5it)] ds. 

Jan Jar 


(3.7) 


If now one recalls the principle of virtual powers according to which a given system is in equilibrium if the power 
of internal forces is equal to the power of external forces, it is straightforward that the expression (3.6) naturally 
suggests which is the correct expression for the power of external forces that a second gradient continuum may 
sustain, namely: 


V e 


= f /® xt Suj dv + f f® xt Suj da + f < 7 ® xt Su™ da + f [7r| xt 5itj] ds, (3-8) 

J n Jan Jan Jar 
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where / ext are external bulk forces (expending power on displacement), f ext are external surface forces (expend¬ 
ing power on displacement), g ext are external surface double-forces (expending power on the normal derivative 
of displacement) and 7 r ext are external line forces (expending power on displacement). Imposing that 

pint + pext = q (3.9) 

and localizing, one can get the strong form of the equations of motion and associated boundary conditions for 
a second gradient continuum. 

Therefore, the equilibrium equation for a second gradient continuum is 

Div (a — Div m) + / ext = 0. (3.10) 

This set of partial differential equation can be complemented with the following boundary conditions: 

• Strongly independent, second gradient, geometric boundary conditions 

u | r = zt ext , 

Vu-n| r = a ext , 

for given functions tt ext , a ext : R 3 —> R 3 on the portion T of the boundary. 

• Strongly independent, second gradient, traction boundary conditions 

Traction boundary condition on dfl \ T: 

(er — Div m) • n — V [(m ■ n) ■ T] : = t ext , traction (3 be) 

(m • n) • «| an yp = g ext , “double force traction” (3 be) (3.12) 

for prescribed functions f ext , g ext : R 3 —> R 3 at the boundary. 

Traction boundary condition on 9T: 

[(m • n) • ^] | ar = 7 r ext “line force traction” (3 be) (3.13) 

for a prescribed function 7 r ext : R 3 —>■ R 3 on <9T. 

We want to stress the fact that in the framework of a second gradient theory the test functions that can 
be arbitrarily assigned on the boundary dfl are the virtual displacement Su and the normal derivative of the 
virtual displacement V (Su) ■ n. This means that one has 3 + 3 = 6 independent geometric boundary conditions 
that can be assigned on the boundary of the considered second gradient medium. Analogously one can think 
to assign 3 + 3 = 6 traction conditions on the force (in duality of Su) and double force (in duality of V(<5u) • n) 
respectively. Hence, in a complete second gradient theory 6 independent scalar conditions must be assigned on 
the boundary in order to have a well-posed problem. 

3.2 The indeterminate couple stress model viewed as a subclass of the second 
gradient elasticity model 

As in the previous case, we consider a particular strain energy density of the type 

W = Wii n (Vu) + Wc Urv (V[axl(skewVu)]) = Wn n (Vw) + W C urv(Vcurlu), 

where Wi; n (Vu) is given in eq. ( |2.2| ), while the curvature energy W CU rv(Vcurlu) also discussed in the previous 
section, is given by 


(3bc) 

(3bc) 


(3.11) 


Wcurv(Vcurlu) 


= “J" II sym Vcurlu || 2 + -^ || skew Vcurl u || 2 

=: y II S II 2 + ^WM \ 2 = Y SlmSlm + 


(3.14) 
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where we set 


S pq := (sym Vcurl it) p9 = €prsUs ’ rq + £qrsUs ' rp , 


A pq := (skew Vcurl u) pq = Cprt,Ua ' rg CgraUfl ’ rp . 


(3.15) 

This decomposition of the curvature energy is equivalent to that which can be found in Mindlin and Tiersten 


and presented by us in eq. ( 1 . 81 . 


Regarding the curvature energy (3.141 as a particular case of second gradient energy, we can directly calculate 

(3.16) 


the particular form of the third order hyperstress tensor as 

cWcurv 


tTy/i' — 


1^11'curv dS pq ( 0 ilcurv UA pq 0C\ ^ ()S pq | (X ‘2 , 3 A pq 


3Ui^j k 3S pq dUijk 3A pq Ov.'jj X 


2 pq duijk 2 pq dui. 


,jk 


It can be checked that, from (3.15), one gets 


Q — o (. e pji^qk + tqjifipk) ■. 

OUijk 2 


dA pq _ 1 x __ ^ 

q u — 2 CqjiOpk) ■ 


Replacing these expressions in (3.16l, using the definitions (3.15) together with the identities 


- S ,r S'j a SncSi‘ 


and 


tpji^krs — d pk Sj r Si s d pk Sj s Si r d pr Sj k 6i s T d pr 6j s 8i k T d ps Sj k Si r d ps Sj r 6i k: 
the fact that S p k = S kp and A pk = —A kp and simplifying gives 


~ Ql , _ Q \ r A , . «1 ^ «2 . 

ffl-ijk — , ytpjibpk I tqji^kq) ' 7“ \£pji.A pk £qjiA kq ) € p jiD pk * ,, €pjiA pk 


OL2 

~4 


Ui 


Ct2 


(3.17) 


_ / , \ _j_ ^2 / \ 

— €pji\€prs'U j s,rk i €krs'U j s,rp) i €pji\€prs'U j s,rk ^krs^s^rp) 

OL\ 1 

= i,jk ^j,ik) H - ^ (^1 ^ 2 ) 'U j i,pp$jk H - ^j,pp^ik 'U j p,jp&ik\ • 

Such particular expression of the third order hyperstress tensor can be also written in compact form as 

m = ai V [skew (S7u)] + ^ (ai — 0 : 2 ) [V (div n)0l- Div (Vu) 0 1] 

+ ^ («i — ct 2 ) (Div (Vzt) ® 1) T - (V (divu) <g) 1) T 

where we denote by the superscript T 12 the transposition over the two first indices of the considered third order 
tensors. With such definition of the third order hyperstress tensor m one can now write the principle of virtual 
powers for the considered particular case in the form 


(3.18) 


/ (<7jj - m ijk ,k) , Sui dv- 
In Jan 


((Jij m rjk[ k ) Tlj (TpjXWjjl; Tip ) q 1 hp 


Sui da 


/ (m.ij k n k nj) Su™ da - [m^*, n k UjSuij ds 

Ian Jar 


(3.19) 


= — f /, ext Su t dv — f tf^dui da — f m| xt 5u" da— f 7r° xt <5uj ds . 
J n Jan Jan Jar 


We have to remark that the term 


_ QU 

(iTl ijk n k fl : j ) dl.L t — (U/jf- Vj ,rf: )a. k a.j 


^ (tr 1 Oi 2 ) (. U p ' /p 7ljTlj Ui^ppUjUj H - Uj^ pp TljTli 'UpjpTljTli^j 


Su 7 


(3.20) 
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is vanishing for some particular choices of the indices. In particular, if, for the sake of simplicity, one 
considers the introduced quantities to be all expressed in the local orthonormal basis (n, t, v}, then the 
aforementioned term can be rewritten as 


(m ijknkrij) Su " = 


a i 


0 (w^, 11 Wl.il) T . (cil Oi 2) (Up : ip Ui^pp Ui^ppTli Up^lpTli) 


su ? 


It can be easily checked that such term is vanishing when i = 1. More precisely, we are saying that the 
normal component of the normal derivative <ht" does not contribute to the power of internal forces when 
considering the indeterminate couple stress model. This is equivalent to say that indeed only 2 geometric 
boundary conditions can be imposed on the normal derivative of virtual displacement or, equivalently, on 
its “traction” counterpart which is the double force. 

Hence, the governing equations of the considered system can also be formally written in the form 


Div(cr — Divm) + f ext = 0, in duality of Su 
together with the following boundary conditions induced by (3.7 


(3.21) 


• Strongly independent, geometric boundary conditions for the couple stress model on T (as 
derived by a full-gradient model) 


u 


r 


T -Vu ■ n 


u ex \ 

T ■ a ext , 


(3.22) 


for given functions u ext ,a ext : R 3 —> M 3 at the boundary. 

• Strongly independent, traction boundary conditions on dfl \ T (as derived by a full-gradient 
model) 


(u — Divm) • n — V [(m • n) ■ T] : T = f ext , in duality of Su (3.23) 

T ■ [(m ■ n) ■ n] = T ■ g ext , in duality of T ■ ( VSu ) ■ n 

for prescribed functions f ext , g ext : R 3 —> R 3 at the boundary. 

Traction boundary condition on the curve dr: 

[(m • n) ■ uj = 7 r ext , (3.24) 

for a prescribed function 7r ext : R 3 — > R 3 on dr. 


3.3 


Reduction from the third order hyperstress tensor m to Mindlin’s second order 
couple stress tensor m 


We want to prove here that the equations (3.101 and the traction boundary conditions (3.22)-(3.23) can be 
equivalently rewritten using Mindlin’s second order couple stress tensor 


~ ai+a2 , «i — a 2 , , , t 

m= --- Vcurirt-|--- (Vcuriu) 

= nr devsym(Vcurlu) + a .2 skew(Vcurlu) 

= 2aq devsym(V axl(skew Vw)) + 2a 2 skew(V axl(skew Vu)), 

^ Oi 1 + Oi-2 Ui — Q!2 

Wli/ = — tijk'U’kJl T tmjkMk.jii 


Mindlin (3.25) 

equivalent form 1 
equivalent form 2 

index format 


3 We recall that in the considered couple stress model expressed in the framework of a full second gradient theory the constitutive 
form for m is given in eq. (|3. 1 7[) or equivalently {|3.18[i. 
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instead of the third order tensor given in eq. (3.181. Such second order hyper-stress tensor has been introduced 
by Mindlin and Tiersten |59] and we have shown in a previous section that it can be obtained by means of a 
direct variational approach that does not need the introduction of the third order couple stress tensor m (see 

eq. @). 

In order to be able to set up such equivalence, we have t o rem ark that, for the choices (3.181 and ( 3.25| ) of 
m and to, the following properties are verified (see Appendix A.2 for detailed calculations^ 


1 ~ l - Ot 2 

Div ^ antiDiv = ---A(skewVu), 

f3x3x3 TO3x3 


(3.26) 


and 


V [(m • n) ■ T] : T = -V[anti ( fh-n)-T]:T , (3.27) 

T ■ [(m ■ n) ■ n] = -T ■ anti(m • n) ■ n, (3.28) 

[(m • n) ■ v\ = ^[[anti(m • n) • n] ■ v\ , (3.29) 


where 


— 1 _ (oti — Q(2) 

m • n = -anti (to ■ n) = a i [V (skew Vu)] • n + - --—— skew [V (Div u) ® n — Div (Vu) ® n]. (3.30) 

Clearly, based upon such relationships, we can recognize the following equivalent forms for the bulk equations 


1 


Div(cr — Divrh) + / ext = 0 <t=> Div — - anti Div fh J + / ext = 0, 

in duality of Su 

together with the following equivalent forms of the traction boundary conditions 


(3.31) 


1 


(a — Divm) • n — V [(m ■ n) ■ T] : T = t ex <t=> — - anti Div to ) ■ n — -V[anti (to ■ n) ■ T] : T = t ex (3.32) 

in duality of Su, 


1 . 


1 . 


T ■ [(m • n) ■ n] =T ■ g ext -T ■ anti(TO • n) ■ n = T ■ g ext 


(3.33) 

in duality of T ■ (VSu) ■ n , 

and finally the following equivalent conditions on the boundary of the boundary <9r where traction is assigned 


(m • n) ■ v\ = 7 r ext <^> ^Janti(TO • n) ■ v\ = 7r ext 
in duality of Su. 


(3.34) 


3.4 A direct way to obtain strongly independent boundary conditions in the 
indeterminate couple model 

Let us consider again the energy 

W = Wiin(Vu) + Wcurv(VcUrl u), 


with Wy ln and Wc Urv defined in equations (2.2) and for which different equivalent forms of the curvature energy 
have been given in eq. (1.8 1 . To the sake of completeness, we derive in this section the equations of motion and 


4 Using a classical notation A = DivV is the Laplacian operator. 
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associated boundary conditions of the couple stress model by directly computing the first variation of the action 
functional associated to the considered energy, without noticing that such energy is indeed a very particular 
case of a second gradient energy. This procedure follows what was done by Mindlin and Tiersten [59] and it 
is presented in Section [2] Here, as done in [59], the curvature energy is regarded as a function of the second 
order tensor Vcurlw, instead that of the third order tensor VVu. The difference of the calculation that we 
present here, with respect to what is done by Mindlin and Tiersten, is that we proceed further in the process 
of integration by parts up to the point of getting strongly independent quantities on the boundary. As it is 
shown in Section |3.3[ the two approaches can be considered to be finally equivalent, provided that a suitable 
identification of the second and third order tensors appearing in the governing equations is performed. 

As usual, the power of internal actions is given by the first variation of the action functional which can be 
directly computed as 


V mt = 5A = —6 


Jn L 
,dW yin 
' dVu 


Wlin(Vtt) + W curv (Vcurl u) 

,dw curv 


, dVu ) dv — 


dW ,; 


lin 


g u Sutj dv - 


in uu i,j 


JnA dS 
f dW cur 

L -Hs- ss « dv - 


dv 
SS) dv — 


,dW r , 


dA 


<5A) dv 


(3.35) 


dW„ 


6An dv. 


■n dA tj - 10 


Using the expression of Wc Urv given in (3.141 and then the definitions (3.151 for S and A together with the 
properties of Levi-Civita symbols, it can be checked that 


dW a 


in dS tj " lJ 
Oq 


SSa dv — 


dW a 


6An dv 


■n dAij "* ZJ 

«2 


=-— / SinSSa dv - — And An dv = 


4 

aq 


in 


in 


j (ti Pq u q , pj + e jpq u q , pi ) (e irs du s , rj + e jrs du s , ri ) dv 

Jn 

ex 2 r 

0 / (^ipqllq.pj ^jpqU’q.pi') ( t'irsdlls.r'j ^jrsdlls,ri ) dv 

n Jn 

(ax + a 2 ) 


£ipq£irsll'q,pj3lls,rj T Hjpq^jrsHq.pi 


, P idVj S ^ r i) dv 


(aq — a 2 ) 

8 

(oil + a 2 ) 


.. 

I [(fipr^qs $ps$qr) 'U , q,pid'U , s,ri\ dv 

Jn 


(ai — a 2 ) 


((dijSp r 5q S dijSp S Sq r dj r SpiSq S ~b Sj r Sp 3 Sqi Sjsdpidqr djs^pr^qi') Uq^piSUs,rj') dv 


(ax + a 2 ) f _ Ur s ) du Sir i dv 

L (<“- - ++ <** - * 


4 

(aq — a 2 ) 


— 2 ^ (Hs.r Hr,s') x dv ^ ~ J Vsyi') x du s , rr A (u r y Uiyr) ^ dll s ^ r Jj dv . 

Recalling also the results for the variation of the classical first gradient term given in (2.41 this last relation 


implies that the internal actions (3.351 can be rewritten as 


V m = dA = - / (/i(uij + u jt i) + Xu k ,kdij) duij dv -- / SijdSij -- / A^dAijdv 

Jn 4 Jyi 4 Jq 

= — (H (uij + Uj ti ) + A u k ,k5ij) duij dv -/ (u SjT . - u r , s ) i du StTi dv 

Jn z Jn 
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(ai — < 22 ) 
4 


'n L 


(^i,s //s,i) ^ hu s ^ rr 4“ (drj ^2,r) ,y hu s ^ r 


dv. 


Suitably integrating by parts we can hence write 
■p mt = <L4 = / ai^rijSui da+ jSut dv 


/an 


a-i , \ , (<*i — ^ 2 ) 

— (u Sjr — Ur,s) ;ii H | 


(5it Sir dv 


1 an 


^1 / \ 1 

\^s,r '^'r,s) j i H - ^ 


^s,i) i r (^r,i ^> r ),is) 

'U'S’i) i Tl r ~b ('U r > ) 2 ^i,r) 1 


(3.36) 


Stic r da 


- / ((Tij - Tij) rijSui da + (<x^- - ) . 5tij dv- B i:j Su id da , 

Jan Jn Jan 


where we set 


(Tij — (// ('U'ijj 4 “ Hj.i) 4 “ Ukj : Sij ) , 

ai l„. „, \ ( a l — a2 ) 

\Uijpp — Uj,ipp) ^ 


( V’i.ppj u j,ppi ) 


(ai + a 2 ) 


( u i,j %,*) • 


4 V~M “JiV.pp: 


and 


o _ ai / v „ . (ai-a 2 ) 

Bij — 2 ( U *>J u j,i),p n P+ 4 


((w p ,i u *,p) , p n j + ( u j,p u p,j), P n *) ■ 


With reference to eqs. ( 3.26[ ) and ( |3.30 >, it can be recognized that 
1 ___ ai -(— c^2 

r = Divm = -anti Div to = ---A(skewVu), 

_ 1 _ (ax — 0 ^ 2 ) 

B = m • n = -anti (m • n) = cty [V (skew Vzt)] • n + --—- skew [V (Div it) Cg) n — Div (Vu) <g> n] 


(3.37) 

(3.38) 

(3.39) 


with m and m given in (3.181 and (3.251 respectively. 


Remark 3.1. We explicitly remark at this point (and we will point it out more precisely in the next section) 
that the results presented by Mindlin and Tiersten m are compatible with a variational procedure which stops 
at this point (eq. (3.36 )) without proceeding further in the process of integration by parts. 

Indeed, in the view of proceeding towards the determination of strongly independent virtual variations, the 
last term in the expression (3.36) of the power of internal forces can still be manipulated according to the 
procedure (3.51 of surface integration by parts, so that one finally gets 

\BijVp8ui\ds - / (T p jBij) h T hp Suida+ / (B t .jnj) Su™ da. 

’ Jan ’ Jan 


/an 


BijSuij da — 


i dr 


Hence, supposing that the virtual displacement is continuous through the curves dr, the power of internal forces 
of the couple stress model calculated by means of a direct approach reads 


= / {(Tij - Tij) j Sui dv- 


I I (C 


z"?a) /f 7 (l Bjj ) y, hfij) Suida j (BijUj) Su i da j | Bj.j iy, Jh//? ds . 


ion 


la r 


It has already been proven in Subsection 3.2 that only the tangent part of the normal derivative 5u™ contributes 
to the power of internal actions when considering the indeterminate couple-stress model, so that the power of 
internal actions can be finally written as 


'P = J ^ ~ Tij ^’ j ^ dV ~ 

/ | B’ij Vp J dzq ds. 

Jar 


/an L 


(<Tij - Tij ) rij - ( TpjBij ) h T h p Sui da- ( T ip Bpjnj ) (T ih 6v%) da 


i an 


(3.40) 
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or equivalently, in compact form 


-pmt _ f (Di v (cr — r) Su ) dv — I ([(<7 — r) • n — (V (B ■ T)) : T] ,5u) da (3-41) 

J a Jan 

— j ((T ■ B ■ n) ,T ■ S(\7u ■ n)) da — f {\B ■ i/J, Su) ds, 

Jdn Jdr 


where we recall once again that the tensors r and B are given by eqs. (3.381, (3.391. 

Considering the power of external actions to take the form (3.8) imposing V + ’P ext = 0 and localizing, one 


gets the bulk equations and associated traction boundary conditions for the couple stress model by means of a 
direct approach 


Div (cr — r) + / ext = 0 in duality of Su , 

together with the following traction boundary conditions on the portion of the boundary dfl \ T 


(3.42) 


(a - t) ■ n — [V (B ■ T)] : T = t e 


in duality of Su , 
in duality of (VSu) ■ n 


T ■ B ■ n = T ■ g ext 

and finally the following condition on the boundary of the boundary 9r where traction is assigned 

\B ■ v\= 7r ext in duality of Su. 


(3.43) 

(3.44) 


(3.45) 


Given the identification of the tensors r and B with the tensors m and m as specified in eqs. ( 3.38| ), (3.39), 
the bulk equations and traction boundary conditions (3.42)-(3.45 I as derived by means of a direct approach are 


completely equivalent to eqs. (3.31)-(3.34). 


3.5 The geometric and traction, strongly independent, boundary conditions for 
the indeterminate couple stress model 

We have proven up to now that, independently of the method that one wants to choose to obtain the correct 
set of bulk equations and associated boundary conditions, passing through a full second gradient approach or a 
direct approach based on second order tensors instead the third order ones, one finally arrives at the following 


complete set of boundary conditions which can be used to complement the bulk equilibrium equation (3.421 of 
the couple stress model: 


3.5.1 Geometric (essential or kinematical), strongly independent, boundary conditions on T 


u = zi ext (3 be) (3.46) 

(1 — n ® n) ■ (Vu • n) = (1 — n <S> n) ■ a ext (2 be) 

where u ext ,a ext : R 3 —> R 3 are prescribed functions on the subportion T of the boundary dfl, where kinematical 
boundary conditions are assigned. 


3.5.2 Traction, strongly independent, boundary conditions on dfl \ T 


Correspondingly to the geometric boundary conditions, we may prescribe the following traction boundary 
conditions based on (3.431 (or equivalently (3.321) 


(<r — t) • n — |V[anti(m • n) • T] : T 

_ j.ext ^ 

’ , l 

on dQ. \ r 

(3 be) 

(1 — n (g) n) ■ anti(fh • n) ■ n 

= (1 — n g) n) ■ g ext , J 

(2 be) 

\ [anti(fh • n) ■ v\ 

= 7T eXt , 

on dr 

(3 be) 


(3.47) 
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where t ext ,g ext : R 3 —> R 3 are prescribed functions on <912 \ T, while 7 r ext is prescribed on <9T and leads to 3 
boundary conditions on the curve <917. 

It can be shown (see Appendix A .6 for the proof of the needed identities (A.31) and (A.32)) that such set 
of traction boundary conditions can be ulteriorly simplified in the following form 


(a — t) ■ n — |V[anti(m • n) • T) : T = t e 


(1 — n ® n) ■ anti((l — n (g> n) ■ m ■ n) ■ n = (1 — n (g) n) ■ g ext , 

\ [anti((l — n ® n) ■ rh ■ n) ■ u\ = 7 r ext , 


on <9f l \ r 


on DV 


(3 be) 
(2 be) 
(3 be) 


(3.48) 


where t ext ,g ext : R 3 —> R 3 are prescribed functions on <9f l \ T, while 7 r ext is prescribed on <9T and leads to 3 


y.ext ~ext 

6 i y 

boundary conditions on the curve <917. 


4 Assessment of the strongly independent boundary conditions for 
the indeterminate couple stress model in a form directly compara¬ 
ble to Mindlin and Tiersten’s ones 


Given that the bulk equations (3.421 that we obtained are the same as Mindlin and Tiersten’s ones, the delicate 
point is now to compare our boundary conditions (3.46)-(3.47) with those provided by Mindlin and Tiersten 
in [59| . If a proof of the equivalence of the purely kinematical boundary conditions (3.461 with those proposed 


by Mindlin and Tiersten has already been provided in Lemma |2.3[ the equivalence between traction boundary 
conditions as derived with our and Mindlin’s approach is not straightforward. This is why we need here to 
rewrite the boundary conditions (3.481 in a suitable form. 


4.1 Towards a direct comparison with Mindlin’s traction boundary conditions 


In order to be able to directly compare the traction boundary conditions for the indeterminate couple stress 
model which we obtained both passing through a second gradient theory and by means of a direct approach 
with those proposed by Mindlin, we need to rewrite our equations in a suitable form. In this section we show 
some calculations which are needed in order to reach this goal. 

Proposition 4 . 1 . For all m £ R 3x3 and for all smooth surfaces S = {(x\,X2, £3) £ K 3 | F(xi, X2, X3) = 0 }, 
F : uj C M 3 — > R 3 of class C 2 , the following identity is satisfied: 

^V[anti (m • n) ■ T] : T = ^n x V [ (n, (sym rh) ■ n )] + [anti (T ■ m • n) ■ T] :T. (4.1) 

Proof. The proof is included in Appendix |A.7| 0 


Remark 4.2. From the above proposition, it follows that the first of the boundary conditions ( 3.47| ) (or equiv¬ 
alently (3.32),) can be finally re-written in the form 


(a — t) • n — x [V (( n , (sym?n) • n))] — [anti (T ■ m ■ n) ■ T] : T = t ext . 


(4.2) 


In Section |2.2| we have recalled the argument of Mindlin and Tiersten and we have remarked, see (2.17), 


that the term — [anti (T ■ m ■ n) ■ T] : T is absent in their formulation since it remains somehow hidden in 
duality of curl(5u) which is not manipulated further in their formulation. 


4.2 Final form of the strongly independent, geometric and traction boundary 
conditions for the indeterminate couple stress model 

Basing ourselves on the previously results obtained in Subsection |3.3| we can now establish which is the set of 
geometric and traction boundary conditions to be used in the indeterminate couple stress model, alternatively 
to those proposed by Mindlin and Tiersten. As we will better point out in the remainder of this section, the 
boundary conditions that we derive by our direct approach are as legitimate as those proposed by Mindlin and 
Tiersten. Nevertheless, if in one case one can equivalently pass from one set of imposed boundary conditions to 
the other one, such equivalence cannot be stated for the case of mixed boundary conditions. 
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4.2.1 Geometric (kinematical, essential) strongly independent boundary conditions for the in¬ 
determinate couple stress model 

As for the geometric boundary conditions, we recall that one can assign on T C dfl the following conditions 

u = u ext (3 be) 

(1 — n ® n) ■ Vu • n = (1 — n ® n) ■ a ext (2 be) (4.3) 

where u ext ,a ext : R 3 —> R 3 are prescribed functions. Such conditions are the geometric boundary conditions 
which are known to be valid in the framework of second gradient theories, with the peculiarity that here only 
the tangent part of the normal derivative of displacement can be assigned here. 

We have already shown that the fact of assigning the tangent part of Vrt • n is indeed equivalent to assigning 
the tangent part of curlu, so that such set of geometric boundary conditions can be seen to be equivalent to 
Mindlin and Tiersten one’s according to Lemma |2.16| 


4.2.2 Traction strongly independent boundary conditions for the indeterminate couple stress 
model 


As far as the traction boundary conditions are concerned, considering the manipulated form (4.2) of equation 


(3.48 11 , the strongly independent boundary conditions (3.481 for the indeterminate couple stress model can be 


finally rewritten as 


[(a — t) ■ n — |n x V[(n, (symm) • n)] 

— |V[(anti[(l — n ® n) ■ fh ■ n]) ■ T] : T 

= t ext , 

l on <9f l \ T 

(3bC> (4.4) 

(2 be) V ; 

(1 — n (g) n) • anti[(l — n ® n) ■ fh ■ n] ■ n 

= (1 — n (g) n) g ext J 

I 

[anti[(l — n ® n) ■ fh ■ n] ■ v\ 

= 7T eXt , 

on dT 

(3 be), 


where t ext ,g ext : R 3 — >• R 3 are prescribed functions on dfl \ T, while 7r ext : R 3 — >• R 3 is prescribed on dT and 
leads to 3 boundary conditions. 

In this section we have deduced the strongly independent traction boundary conditions which are coming in 
a natural way from second gradient elasticity and we have compared them to those presented by Mindlin and 
Tiersten thus showing their apparent disagreement. 


5 Are Mindlin and Tiersten’s weakly independent boundary condi¬ 
tions equivalent to our strongly independent ones? 

Up to this point, we have shown that the boundary conditions derived by Mindlin and Tiersten (59] for the 
indeterminate couple stress model are not directly superposable to those that we obtain by means of a standard 
variational approach in the spirit of second gradient theories. 

Even if these sets of boundary conditions are formally not the same, they both follow from the same strain 
energy density. The only difference that we can point out in the two approaches is related to the process of 
integration by parts which is perfomed on the action functional based upon the considered strain energy density. 
Indeed, Mindlin and Tiersten’s boundary conditions are only "weakly independent", while those obtained by 
means of our direct approach can be considered to be "strongly independent" in the sense established in 
Subsection 11.41 

To the sake of compacteness, we use in the sequel the following notations for the internal tractions and 
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hypertractions respectively as obtained by Mindlin and Tiersten’s and our approach 


t lnt := ( a — - anti(Divm) ) • B--nx V[(n, (symm) • n}], 


t' nt := (cr — ^ anti(Divm)^ • n — x V[(n, (symm) • n)] 

— ^ V [ anti ( (1 — n®n) -fh-n) • (1 — n ® n) ] : (1 — n (g> n), 

<f nt := (1 — n ® n) ■ fh ■ n, 

g lnt := (1 — n ® n) ■ anti[(l — n ® n) ■ m ■ n] • n = anti[(l — n (g> n) ■ m ■ r] ■ n. 


Mindlin-Tiersten’s formulation 
our formulation 
(5.5) 

Mindlin-Tiersten’s formulation 
our formulation 


In the last equality for g' nt we have used the fact that g lnt is indeed the dual of T ■ V(5u) • n which is a 
vector tangent to the boundary, which is equivalent to say that the normal part of g lnt does not intervene in 
the balance equations. 

To the sake of of simplicity, we summarize the two sets of possible geometric and traction boundary conditions 
as obtained by Mindlin and Tiersten and by ourselves in the following summarizing box 


Table 1. Possible sets of boundary conditions in the indeterminate couple stress model. 



Mindlin and Tiersten 

Our approach 

Geometric 

(I) u = u ext , 

(II) T ■ curl u = T ■ a ext 

(III) u = u ext 

(IV) T ■ Vm ■ n = T ■ a ext 

Traction 

(A) Z int =? xt , 

(B) g int = <? ext 

(C) f int = t ext 

(D) g int = g ext 

Mixed BCs 1 

(I) u = u ext , 

(B) g int = g ext 

(III) u = u ext 
(D) 5 int = g ext 

Mixed BCs 2 

(II) T ■ curlu = T ■ a ext , 
(A) t int = t ext 

(IV) T ■ Vm ■ n = T ■ a ext 
(C) f int = t ext 


The problem now arises to establish the equivalence between analogous sets of boundary conditions in the 
two approaches. Since all the presented boundary conditions arise from the same strain energy density, we would 
naively expect a complete equivalence between the two models. We will instead show that, if a direct equivalence 
can be established in some cases, this is not indeed feasible for all possible sets of boundary conditions that may 
be introduced in couple-stress continua. More particularly, we individuate different possible sets of boundary 
conditions that are allowed in couple-stress continua being compatible with the Principle of Virtual Powers as 
settled in Mindlin’s and Tiersten’s and our approach respectively: 

• Fully geometric boundary conditions. The boundary conditions (I) and (II) (Mindlin and Tiersten) 
or (III) and (IV) (our approach) are simultaneously assigned on the same portion of the boundary. 

• Fully traction boundary conditions. The boundary conditions (A) and (B) (Mindlin and Tiersten) 
or (C) and (D) (our approach) are simultaneously assigned on the same portion of the boundary. 

• Mixed 1: displacement/double-force boundary conditions. The boundary conditions (I) and (B) 

(Mindlin and Tiersten) or (III) and (D) (our approach) are simultaneously assigned on the same portion 
of the boundary. 

• Mixed 2: force//.) 1 (it) boundary conditions. The boundary conditions (II) and (A) (Mindlin and 
Tiersten) or (IV) and (C) (our approach) are simultaneously assigned on the same portion of the boundary. 
We recall that by D x (m) we compactly indicate the operator T ■ curl u when we consider Mindlin and 
Tiersten’s approach or the operator T • Vtt • n when considering our approach. 
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We explicitly remark that, in order to be consistent with the introduced Principle of Virtual Powers, when the 
first sets of conditions is applied on a portion P of the boundary, the second ones must be assigned on the 
portion <9fl\T. Analogously, when assigning the third set of boundary conditions on T, the fourth one must be 
assigned on \ T. 

In the following subsections we carefully study the four introduced cases by establishing whether Mindlin 
and Tiersten’s approach is equivalent with our formulation of the indeterminate couple-stress model. 

5.1 Fully kinematical boundary conditions 

We have already shown (see Lemma |2.3[ ) that it is equivalent to simultaneously assign the displacement and 
the tangential part of its curl or the displacement and the tangential part of its normal derivative on the 
same portion of the boundary. This means that Mindlin and Tiersten’s boundary conditions (I)+(II) are 
completely equivalent to our conditions (III)+(IV) in the sense that one system of equations can be directly 
obtained from the other and vice-versa. We would like to thank an unknown reviewer for pointing out to us 
this equivalency. 

5.2 Fully traction boundary conditions 

We consider here the case in which forces and double forces are simultaneously applied on the same portion 
of the boundary. More particularly, this means that we are simultaneously applying on the same portion of 
the boundary conditions (A) and (B) (Mindlin and Tiersten’s) or conditions (C) and (D) (our approach). 

We start by showing that conditions (B) and (D) are equivalent. To do so, we notice that, given two vectors 
v and n, one can check that, according to definitions |1.3| the following equalities hold 

(antiu)ij rij = —Cijk i’k rij = ~(n x n)i = (v x n)i. 


This means that we can write: 


g mt = anti[(l — n ® n) ■ m ■ n] ■ n = ((1 — n<8> n) ■ m ■ n) x n, 

so that the boundary condition (D) can be rewritten as ((1 — n ® n) • m • n) x n = g ext . Comparing this 
last equation with equation (B), we can finally conclude by direct inspection that equations (B) and (D) are 
equivalent when setting 

g ext = g ext x n. (5.6) 

On the other hand, eq. (C) can be rewritten as 

^er — ^ anti(Divm)^ 'U-^nx V[(n, (symm) • n)] = f ext + [ anti (T ■ m ■ n) ■ T] : T, 

which, considering eq. (B) can also be rewritten as 

(a — ^ anti(Divm)^ • n — x V[(n, (symm) • n)] = f ext + [ anti ( g cxt ) • T] : T 

already known since only the tangetial derivatives are considered 


It is easy to check that this last equation is equivalent 








text = t ext - 1V [ anti (if xt ) • T ] : T. 


(5.7) 


We have thus proved that, in the case of fully traction boundary conditions, given a couple of tractions (f ext , g ext ) 
in Mindlin and Tiersten’s model, one can always “a priori” find a corresponding pair of tractions (t ext ,g ext ) in 
our model such that the two sets of boundary conditions (A)+(B) and (C)+(D) are equivalent thanks to the 
relationships (5.6 1 and (5.71. The converse is clearly also true. 
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5.3 Mixed 1: displacement/double-force boundary conditions 

We treat here the case in which we simultaneously assign the displacement and the double force on the same 
portion of the boundary. This is equivalent to say that one is assigning eqs. (I) and (B) in the Mindlin 
and Tiersten’s approach or (III) and (D) when considering our approach. As already proven in the previous 


subsection, the equivalence between equations (B) and (D) can be obtained by setting the relationship (5.6 1 
between g ext and g ext . Moreover eqs. (I) and (III) are clearly equivalent when u ext = u ext . 


5.4 Mixed 2: fore e/D l (u) boundary conditions 

We have shown up to now that for the three preceeding cases of boundary conditions an a priori equivalence can 
be established between Mindlin and Tiersten’s and our couple stress model. More particularly, this means that, 
given a boundary value problem stemming from Mindlin and Tiersten’s model, we can set up (thanks to suitable 
identifications between tractions and double tractions in the two models) another boundary value problem which 
give rise to the same solution. We show here that the establishment of such an a priori equivalence is not possible 
in this last "Mixed 2" case in which forces and higher derivatives of the displacement field are simultaneously 
assigned on the same portion of the boundary. We start by noticing that if the displacement held is not assigned 


on the boundary, then Lemma 2.31 is not valid any more and, as a consequence, equations (II) and (IV) are no 
longer equivalent. Indeed, starting from the equations (1.211, it is possible to easily deduce that the tangent 


part of the normal derivative and of the curl of the displacement held are respectively given by 

rri t -7 f &ii T du v \ . , f du n duv du T 0u n ^ 

T-Vu-n =\-—, -—, 0 and T-curlu= ---—, ---—, 0 

\OX n OX n ) XOXv ox n ox n OX T 


(5.8) 


When, as in this case, the displacement is not assigned on the boundary, both its normal and tangential 
derivatives are free, so that we cannot establish an a priori equivalence between equations (II) and (IV). The 
same is true when one wants to compare equations (A) and (C). In fact, we can always recognize that eq. (C) 
can be rewritten as: 

(a — - anti(Divm)^ • n — -n x V[(n, (symro) • n)] = t ext + -V [ anti ( g ext ) • T] : T. 


Nevertheless, contrary to the case treated in Subsection 5.2 this last equation cannot be claimed to be equivalent 
to eq. (A) just by setting 


i ext = t ext - ^\7 [anti(T ■ m ■ n) ■ T] : T = f ext — [ anti ( g int (u) ) -T] : T. 


(5.9) 


As a matter of fact, the term g lnt (u ) appearing in equation (5.91 and defined in eq. (5.5 (3 depends on the 


displacement held u through the couple-stress tensor m defined in (3.251. This means that, actually, formula 
(5.9) does not allow to calculate a priori the force t ext in our model which is equivalent to an assigned f ext 


in Mindlin and Tiersten’s model. Of course, one could think to assign, e.g., the boundary conditions (II) and 
(A) on the same boundary, solve the associated boundary value problem so finding its solution u* and then 
calculate the quantities a ext and f ext to be assigned in our model in order to give rise to the same solution it*. 
The converse operation of imposing our boundary conditions (IV) and (C) in our model and then calculate 
a posteriori the quantities t ext and a ext to be assigned in Mindlin and Tiersten’s model to obtain the same 
solution, can also be envisaged. In other words, we are saying that only an a posteriori equivalence is eventually 
possible in the case that the considered force/Z3 1 (it) boundary conditions are applied on the boundary. This 
fact poses, at least, some philosophical problems by giving rise to the question: how to chose among Mindlin 
and Tiersten’s or our mixed boundary conditions?. This indeterminacy leaves many questions open concerning 
the physical transparency of higher gradient theories. 


6 Conclusions 

The present paper gives a comprehensive analysis of the indeterminate couple stress model and of the boundary 
conditions arising in this theory. We have seen the indeterminate couple stress model as a special case of the 
full strain gradient elasticity, we have directly derived the equilibrium equations and the boundary conditions 
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in the same spirit as in the general strain gradient elasticity approach and we have compared our approach with 
that proposed by Mindlin and Tiersten [55] , 

As a balance for the present paper we can state that an apparent inconsistency is found between the classical 
(Mindlin and Tiersten’s) approach and our direct approach to the indeterminate couple stress model. 

Indeed, if an "a priori" equivalence can be found in most cases between the two models, we point out that 
this is not the case when considering "mixed" boundary conditions for which "forces" and suitable combinations 
of first order derivatives of displacement are simultaneously assigned on the same portion of the boundary. 

It turns out that for such particular mixed boundary conditions an "a priori" equivalence cannot be estab¬ 
lished between the two models. 

This fact poses serious conceptual and philosophical problems concerning the transparency of the physical 
meaning of the boundary conditions in higher gradient models. The question remains open whether the use of 
one model would be preferable to the other, at least for the quoted case of mixed boundary conditions. 
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Appendix 


A.l First variation of a second gradient action functional and principle of virtual 
power in compact form 

In this section we basically propose again the calculations concerning the first variation of a second gradient action functional by 
means of a compact notation instead of using Levi-Civita index notation as done instead in Section |3.1[ To this purpose, let us 
consider the second gradient energy ff(Vu, VVu) = W(Vu) + Wcurv(VVu) and the associated action functional in the static case 
(no inertia considered here) 

A = - [ W(Vu, VVu) dv = — [ [W(Vu) + W CUTV (VVu)]dv. 

J n Jn 

The first variation of the action functional can be interpreted as the power of internal actions P int of the considered system and 
can be computed as follows 

V int = 5A = - f (( D Vu W(S7u), VSu) + (D VVu W C urv(Wu), VVSu)) dv, 

Jfl 

Integrating a first time by parts and using the divergence theorem we get 

SA= I (Di v[D^/ u W(Vu)],Su) dv — I (D^/ u W(Vu) ■ n, Su) da 

Jn Jdfi 

— / (D\ 7 ^/ u W C urv (VV u) • n, V Su) da / Div[D^ Vu W C urv(VVu)], VSu) dv. 

JdQ Jn 

Integrating again by parts the last bulk term we get 


SA = — f ({D^ u W(Vu) — Div[DvvuW4urv(VVu)]} • n, 8u) da 
JdQ 

+ / (Div {D Vu W(Vu) - Div[H V VuWcurv(VVu)]} ,5u) dv - / (P V v«^c U rv(VVu) • n, VSu) da, 

Jn J an 

which can also be rewritten as 


SA = - 


/ ({cr — Div[m]} • n, 5u) da + / (Div {cr — Div[m]} , dw) du — / (m • n, VSu) da 
Jan Jn Jan 


(A.l) 


if one sets 

a = D^ u W(Vu), m = Dvvu^curv(VVu). 

Using the brief digression concerning differential geometry (see Sect ion [l. 3 [ ) and recalling the properties we can now ulteriorly 

manipulate the last term in eq. ( |A.1[ ) as follows 

(m • n, \75u) = (m • n, X7Su 1) = (m • n, S78u (T + Q )) = ((m • n) ■ T, X75u) + (m • n, (VSu) • Q) (A.2) 

= ((m • n) • T • T, VSu) + (m • n, (V<5u) • Q) = ( T, T • (m • n) T • VSu) + (m • n, (V5u) • (n ® n)) 

= (T, T • (m • n) T • VSu) + ({n (g> [(V<5u) • n]}m • n, 1) = (T, T • (m • n) T • VSu) + (n 0 {(m • n) T • (V<5u) • n}, 1) 

= (T, T • (m • n) T • VSu) + (n, (m • n) T • (VSu) • n) = ( T, T • (m • n) T • VSu) + ((m • n) • n, (Vdw) • n). 
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We can hence recognize in the last term of this formula that the normal derivative — [(V<5it) • n]^ = Sui^h^h of the 

displacement field appears. As for the other term, it can be still manipulated, suitably integrating by parts for F an open subset 
of dfl and then using the surface divergence theorem ( |1.17| >, so that we can finally write 


I (m ■ n, X78u) da = I (T, V[T • (m • n) T ■ 8u] — V[T • (m • n) T ] • 5u) da + / ((m • n) • n, (V<5u) • n) 

JdQ JdQ JdQ 

= I J((tn • n) T • 8u, „)] ds — / (T, V[T • (m • n) T ] • 5u) da + I ((m ■ n) • n, (V<5u) • n) da 

Jd r JdQ JdQ 

= / [((m • n) • v, <5u)] ds — I ( V[(m • n) • T] : T, 8u) da + / ((m • n) • n, (V8u) • n) da. 

Jd r JdQ JdQ 

The final variation of the second gradient action functional given in (|A.1[), can hence be finally written as 


da 


(A.3) 


5A = / (Div {a — Div[m]} , 5u) dv — I ({cr — Div[m]} • n — V[(m • n) • T] : T, 5u) 

Jq JdQ 


da 


f <C 

JdQ 


m • n) ■ n, ( X78u ) • n) da 


f 1 

J dr 


[((m- n) 


■ v, (5u)] ds, 


(A.4) 


which is directly comparable with eq. GED obtained by means of calculations via Levi-Civita index notation. If now one recalls the 
principle of virtual powers according to which a given system is in equi libriu m if and only if the power of internal forces is equal to 
the power of external forces, it is straightforward that the expression ( |A.4| ) naturally suggests which is the correct expression for 
the power of external forces that a second gradient continuum may sustain, namely: 


^rjext _ 


f (r 

Jq 


*, 8u) dv - 


L 


dQ\T 


(£ ext , Su) da + 


/. 


dQ\r 


(M ext , (V8u) ■ n) da + / (n ext , 8u) ds, 


/ ( 

J dr 


where 

• / ext are external bulk forces (expending power on displacement), 

• £ ext are external surface forces (expending power on displacement), 

• M ext are external surface double-forces (expending power on the normal derivative of displacement) and 

• 7r ext are external line forces (expending power on displacement). 

Imposing that ^p int = —p ext and localizing, one can get the strong form of the equations of motion and associated boundary 
conditions for a second gradient continuum. 


A.2 Some useful relationships between the third order hyperstress tensor and the 
second order couple stress tensor for the indeterminate couple stress model 


Let us consider the third order couple stress tensor 


m = cti V [skew (Vu)] + — (cti — 02 ) [V (Div u) <g> 1 — Div (Vu) <g> 1] 
+ I (ai - an) [(Div(Vu) ® l) Tl2 - (V (Divu) ® l) Tl2 ] , 


together with the second order couple stress tensor 


V curl u - 


2 2 
It can be checked that, suitably deriving eq. (3.171 one gets 

(Divm)ij = Xtlijk,k = ( u i,jkk 'd'jjikk) T ~ (cti Ct2) [Up 
or equivalently in compact form 


V(curlu) T = a 1 devsym(Vcurlu) + 0:2 skew(Vcurlu). 


Cti + OL2 


'iPPJ ■ U'JiPP' 1 u 'P,3P' t l 


( u i,j u j,i) 


,kk ’ 


Divm = 1 —-A(skewVu). 


On the other hand, one has that 


1 r _^ ~ _ 1 ^ ^ ai + Q!2 _ , Oil — OL2 

— [anti DIV 771 J — — — €-ijpTripm,m — ~ ^ijp l ^ e pqk u k,qmm T “ € -mqk' Uj k,qpm 


Cti + OL2 

4 

Cti — Ct2 


(A.5) 


(A-6) 


(A.7) 


(A-8) 


4 

Cti + 

4 

Cti + Ct 2 


(SqjSik Sqj8jk ) Uk,qmm 

( &mi&qjSfcp SmiSqpSfcj -\- SmjSqiSfcp SynjSqpS} 8 rn p 8 qj i 8 ] i j -f- S^pSqjSf^^j '^‘k,qpm 

’P,jpt + 'd'j,ppi H - W p,ipj d,i,ppj 'd'jiipp + Uijpp) 


(\ | Cti Ct2 ^ 


(V/j n Un j ) 
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This relationship, together with ( |A.7|| a nd )A.8| implies the relationship ( |3.26| ). 

As for proving the equalities ( |3.26[ >-( [3.29[ ). we start remarking that 

(m • = (ui,jk n k H - ~ (^1 ^ 2 ) \y>p,ip&jk 'U j i,ppfijk H - ^j,pp^ik ^p,jp^ik\ ^k 

Oq 1 

— ~ \V'i,j ) k ^ (.^1 ^2 ) i^ j p,ip'^ / j Ui^ppl^j ~\~ Uj^ppTl{ Up^jpTli\ , 

which in compact form equivalently reads 


m • n = Qi [V (skew Vu)] • n + — (aq — < 22 ) [V (Div u) 0 n — Div (Vw) 0 n] 

1 'T' 

= — - (aq — Q 2 ) [(V (Divu) 0 n) — (Div (Vw) 0 n)\ 

(Oil — Oi 2 l 

= ai [V (skew Vu)] • n -\ -- skew [V (Div u) 0 n — Div (Vw) 0 n] . 


On the other hand, using the notations introduced in 0 , one has that 


— anti (m • n) 
2 v ' 




Oq + OL 2 


€kpq'U , q,pm ' 


C*q — Q!2 


€mpq'U'q,pk I 


Oil + OL2 

4 

aq — 0:2 
4 

C*q + Q2 


2 2 
(SipSjq SiqSjp')V'q,pm'ft'm 

( ^mi^jp^kq H - ^mi^qj^kp ~\~ ^ip^jm^kq &ip$qj&km ^kp^qi^mj ^mk^pj^qi) ^q^pk^r) 


( aq - a 2 ( . . x 


_ / \ 

— ( U i,j ~ u j,i) 


. Oil ~ OL2 ( x 

m^m T \'U'q,iqTlj 'U'i^ppTlj'lLj^ppTli Uq^jqTli) . 


This relation, when compared to (IA.9I) finally allows to prove that 


_ 1 (oq — a 2 ) 

m • n = -anti (m • n) =ai [V (skew Vu)] • n H---- skew [V (Div it) 0 n — Div (Vu) 0 n] . 

It is clear that, since also in this case m • n = ^anti (m • n), then (|3.26|)-([3.29|) are straightforwardly verified. 


A.3 Some alternative calculations useful to rewrite the governing equations and 
boundary conditions in a form which is directly comparable to Mindlin’s one 

In this subsection we just report some alternative calculations to obtain the same results as before, so that they are not 
uncontournable to the understanding of the main results of the paper. Indeed, the result of the first part of this sections 
can also be re-obtained remarking that 

V[(anti[m • n\) ■ T] — V[(anti[(l — n 0 n) ■ fh ■ n[) ■ T] = V[(anti[(n 0 n) ■ rh • n]) ■ T] 

= V[(anti[n (n, (symm) • n)]) ■ T] = V[ (n, (symm) • n) anti(n) ■ T] (A.9) 

and 

anti(n) (n 0 n) = (anti(n) n) 0 n = (n x n) 0 n = 0, 

( n ® n) anti(n) = — n ® (anti(n) n) = — n 0 (n x n) = 0. (A. 10) 

and also that 

(V[(n,(symm) ■ n) anti(n) • T] : T)i = [( n , (symm) • n) anti(n) (1 — n 0 n) ]ij : k Tjk 
={{n, (symm) • n) anti(n)]<j,*r Tjk — [(n, (symm) ■ n) anti(n) (n 0 n) ]ij : k Tjk 

= [(n, (symm) ■ n) anti(n)]ij, fc T jk (A.ll) 

=[(n, (symm) ■ n)] ]fc [anti(n)]y T jk + (n, (symm) • n) [anti(n)]y,fc T jk 

=[(n, (symm) ■ n)] ]fc [anti(n)]y 5 jk - [(n, (symm) • n)] ik [anti(n)]y rijn k + (n, (symm) • n) [anti(n)]y, fc T jk 
= [(n, (symm) • n)] tk [anti(n)] ife - [{n, (symm) • n)] ik [anti(n)]y njti k + {n, (symm) • n) [anti(n)]y, fe T jk 

= — {n X V[(n, (symm) ■ n)]}i — [(n, (symm) ■ n)],k [anti(n)]y (n 0 n)j k + (n, (symm) ■ n) [anti(n)]y,fc Tj k 

= — {n x V[(n, (symm) ■ n)]}i — [(n, (symm) ■ n)],fc [anti(n) (n 0 n)]ifc + (n, (symm) ■ n) [anti {n)]ij }k T jk 

= - {n x V[(n, (symm) ■ n)]}i + (n, (symm) • n) [anti(n)]y,fc T jk . 


33 



















A.4 The missing steps in Mindlin and Tiersten’s classical approach 

In this section, we present once again the reasoning followed by Mindlin and Tiersten to obtain their set of bulk equations and 
boundary conditions trying to highlight the points in which their approach had to be further developed. 

We start our analysis, remarking that the quantity ( fh-n , (1 — n(gm)• [axl(skew V(5u)]) still does contain contributions performing 
work against Su alone (even though there is a projection 1 — n<g> n involved), which can be assigned arbitrarily and are therefore 
somehow related to independent variation Su. This case is not similar with the Cosserat theory in which we assume a priori that 
displacement u and microrotation A E so(3) are independent kinematical degrees of freedom. On the other hand, the indeterminate 
couple stress model is not simply obtained as a constraint Cosserat model m, i-e. assuming that A = axlskewVu. In the 
indeterminate couple stress model the only independent kinematical degree of freedom is u. We believe that the indeterminate 
couple stress model constructed as a constraint Cosserat model represents only an approximation of the indeterminate couple stress 
model, in the sense that the boundary conditions are not correctly and completely considered. 

Indeed, using the projectors (T = 1 — n 0 n and Q = n<g> n) we obtain 

([(1 — n <g> n) • fh • n], curl Su) = (anti[(l — n (g> n) • fh • n] ■ 1, \7Su) = (anti[(l — n 0 n) • fh • n] • (T + Q), X78u) 

= (anti[(l — n <g> n) • fh • n] ■ T, X78u) + (anti[(l — n 0 n) • rh • n] ■ Q, V<5u) 

= (anti[(l — n 0 ri) • rh • n\ ■ T • T, \78u) + (anti[(l — n 0 n) • m • n]V<5u • Q) da (A.12) 

= (T, T • anti[(l — n (g) n) ■ rh • n] T • \78u) + (anti[(l — n ® n) • rh ■ n], S/8u • Q) 

= — (T, T ■ anti[(l — n 0 n) ■ rh ■ n\ • V5u) + (anti[(l — n (g> n) • rh • n], V5u • Q). 

On the other hand we have 


and, moreover, 


Hence 


(anti[(l — n <g> n) - fh • n], \75u • Q) = (anti[(l — n (g> n) • rh • n], X75u • (n <g> n )) 

= (anti[(l — n (g> n) • fh • n] ■ n, S7Su • n) 

= ((1 — n<g> n) ■ anti[(l — n <g> n) • rh • n] • n, \/5u • n) 
+ ((n <S> ri) • anti[(l — n <g> n) ■ fh - n\ • n, WSu ■ n) 

= ((1 — n<S> n) ■ anti[(l — n(& ri) • fh - n\ ■ n,\75u • n) 

— (n <S> [anti[(l — n (g) ri) • fh ■ n] • n] ■ n, V 8u ■ n) 

= ((1 — n <S> ri) ■ anti[(l — n® ri) • fh • n\ ■ n,\75u ■ n) 

— (n([anti(m ■ n) ■ n], n), V<5u • ri) 

= ((1 — n <S> ri) ■ anti[(l — n® ri) • fh • n\ ■ n,\75u • n) 

— (anti[(l — n <g> n) ■ fh • n] • n, n) (n, \75u ■ ri) 


(anti[(l — n 0 n) • m • n\ -n, n) = 0. 

€so( 3) 


(anti[(l — n (g) n) • m • n\, VSu • Q) = ((1 — n<g> ri) • anti[(l — n <g> n) ■ m • n\ • n, V<5u • n). 


Thus, we deduce 


(m • 7i, (1 — n(g) n) • curl 5u ) = — ( T, T • (anti[(l — n (g> ri) • m • n]) • V5u) 


not completely independent 
second order variation 


((1 — n (g) n) • anti[(l — n <g> n) • m • n] • n, X75u • n ). 


completely Su -independent second 
order normal variation of gradient 


Since 


{V[T • (anti[(l — n<g> ri) ■ m ■ n]) • = {T • (anti[(l — n<g> n) • m • n]) • 8u}i ^ 

= {(T- (anti[(l - n®n) • fh ■ n])) ^{Srijj}^ 

= (T • (anti[(l — n <g n) ■ m • + ( T ' (antiK 1 ~ n 0 n) • m ■ 

= (T ■ (anti[(l — n®n) ■ fh • n ])) ij k(^ u )j + ■ (anti[(l — n®n) ■ fh ■ 

= (T • (anti[(l — n®n) ■ fh • n ])) ij • (anti[(l — nf&n) ■ rh ■ n]) • V<5u}i*. 


(A.13) 


(A.14) 


(A.15) 

(A.16) 


(A.17) 
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using the surface divergence theorenj^Jwe obtain 

/ (fh • n, (1 — n (g) n) • curl Su) da = — / (T, V[T • (anti[(l — n 0 n) • m • n]) • Su]) dc 

JdQ. JdQ 

+ [ T ik (T • (anti[(l - n (g) n) • m • n])) k (Su ) 
Van 


da 


surface 

divergence 


+ / ((1 — n (g) n) • anti[(l — n (g) n) • m ■ n\ • n, \78u ■ n) da 

Jan 

— / [(anti[(l — n (g) n) • m • n])Su , i/)] ds 

Jar 

+ / (T • (anti[(l — n<g> n) • m • n]))^. ,fc Tik (8u)j da 
Jan 

+ / ((1 — n <S> n) • anti[(l — n<g> n) • fh • n] • n, V<5w • n) da 

Jan 

= / ((1 — n (g) n) • anti[(l - n 0 n) • m • n] ■ n,'WSu • n) da 

Jan 

+ / (T • (anti[(l — n<g> n) • rh • n])) i . k (du)j da 

Jan XJ ' 

+ / ([anti[(l — n 0 n) • m • n] • i/], <5u) ds. 

Jar 

On the other hand, we have 

(T • (anti[(l — n (g) n) • fh • n]))^. = T^(anti[(l — n 0 n) • fh • n])^- 

= — T^(anti[(l — n (g) n) • m • n])jz = —T^(anti[(l — n<g> n) • fh • n])ji 
= —(anti[(l — n<S> n) • rh • = —{(anti[(l — n<g> n) - rh • n]) • T}ji. 

Hence, we deduce 

(T • (anti[(l - n (g) n) • fh • n ])). jk T ik = -{(anti[(l - ra <g> n) • m • ra]) • T ik 

= — (V[(anti[(l — n<S> n) • rh ■ n]) • T] : T) .. 

In view of the above computations, we deduce 

— / ((cr — t) ■ n, Su) da — / (fh • n, axl(skew \7Su)) da 

Jan Jan 

= — f ( (cr — r) • n - n x V[(n, (symm) • n)] ,5u)da — f (fh ■ n, (1 — n 0 n) ■ [axl(skew Vdu)]) da 

Jan v 2 Jan '-v--' 


(A.18) 


(A.19) 


(A.20) 


__ /an 

classical first boundary term of Mindlin, Gao, Yang, etc. c * £ 

...must be split further to obtain strongly independent variations... 


isical second boundary term of 
Mindlin, Gao, Yang, etc. 


(A.21) 


r i _ i _ i _ 

/ ((cr-antiDiv[m]) • n - n x V[(n, (symm) • n)]-V[(anti[(l — n<g> n) • m • n]) • T] : T, Su) 

JdCi 2 2 2 

— f ((1 — n 0 n) • anti[(l — n <g> n) • fh • n\ • n, V (5-u • n) da 
2 s --^- / 


da 


completely 6independent second 
order normal variation of gradient 


— f ([anti[(l — n ® n) • m • n] • z/], <5u) ds 

2 Jar 


(A.22) 


for all variations Su E C°°(r2). Note that <9(<9f2 \ T) = dr. Hence, there are indeed two terms 

(a — — antiDiv[m]) • n — —n x V[(n, (symm) • n)] — — V[(anti[(l — n (g) n) ■ fh • n]) • T] : T 

and 

[anti[(l — n <S> n) • m • n] ■ v\ 
which perform work against Su, while only the term 

(1 — n ® n) • anti[(l — n <g> n) • fh • n] ■ n 

is related solely to the independent second order normal variation of the gradient \/Su • n. This split of the boundary condition is 
not the one as obtained e.g. by Gao and Park ED among others and seems to be entirely new in the context of the couple stress 
model. 


5 The surface divergence theorem in this context reads /^(V(T u), T) da = ^)] ds , see (1.151. 
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A.5 The proof of Lemma 2.3 


In other words, this lemma implies that, for any tangential vector field on f C 


(curlu, Ti) | r = known J ( (Vu-n,Tj' 

Since u E C°°(Q), on the one hand we obtain 

u| r = u ext ^ ^ also (curl u, n )| r = known), 


= known, 


(curl u,Tj) 


= known 


curl u |p = known 
u\ r = u ext , 
skewVu | r = known. 


On the other hand we deduce 


u|p = u ext (=>- also V u. Ti |p = known) 


(Vu ■ n,Ti)| r = (b 0 ,Ti) (n,(Vu) T ■ Ti)\ r = (b 0 ,Ti) 


U „ = U 


known 

known, 

known, 

known 


(Vu) T -Ti| r = known, 
(Vu) • Tj | p = known 
Until now, we have obtained 


(skewVu) • Ti\ r 
(sym V u) ■ n |p 


_ .,ext 


(1 — n (g) n) ■ curl u |p = known 


— -,ext 


v« ■ Tj|p 
(skew Vn) ■ T»| r 
(skew Vu) ■ n | r 


u 

known, 

known, 

known 


(Vli.Ti,Ti)|p 
(Ti,(Vn) T .Ti>| r 

(V u.Ti,n) |p 
(x„(Vu) r -Ti) |p 

= u ext , 

= known, 

= known. 


u| p = u ext , 
Vu • Tj | p = known, 
skew Vu | r = known 


u\ r = u ext , 


Vu • 


(sym V u) ■ n | 


= known, 
= known, 


while 


(1 — n (g) n) • Vu • nL = known 


(skew Vu) • Ti| r = known, 
(skewVu)-nL = known, 


“Ian — ’ 

Vu • Tj |p = known, 


(skewVu) • Ti 
(sym Vu) • ^ 


= known, 
= known. 


The only one boundary condition which appears in the last term of ( [A 
(skewVu) • n| r = known. However, (skewVu) • Tj| r = known implies also 

((skew Vu) • n, n) = 0 always, 

((skewVu) • n, Ti) = —(n, (skewVu) • Tj) = known. 
Therefore, by eliminating the redundant information we obtain 

u| p = u ext , 


,26| ) but does not appear in the last term 
(skew Vu) • n| r = known, since 


(skewVu) • n = known. 


(l-n®n)' curluL = known 


_ „,ext 


and 


u |p = u ext , 


(1 — n (g) n) • Vu • n| r = known 
and the proof is complete. 


(skew Vu) • Ti I 


(skewVu) • Tj | 


u 

known, 


_ „ f ext 


u 

known, 


(A.23) 

(A.24) 

(A.25) 


(A.26) 

(A.27) 

of i fOT) is 
(A.28) 

(A.29) 

(A.30) 
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A.6 The proof of some identities to ulteriorly simplify the traction boundary 
conditions 

In this appendix we prove some identities which must be used in order to show the equivalence between some traction boundary 
conditions. We start by checking that the following identity holds 


— (1 — n 0 n) • anti[(l — n 0 n) ■ m ■ n] ■ n = — (1 — n 0 n) • anti(m • n) • n. 


(A.31) 

z 

In other words, we have to check if 

(1 — n 0 n) • anti[(n 0 n) • rh ■ n] ■ n = 0. 

But this fact follows immediately from 

anti[(n 0 n) • rh ■ n] • n = anti [n(n, rh • n)] • n = (n, rh • n) anti[n] • n = (n, rh ■ n) n x n = (n, rh ■ n) 0 = 0. 

The final step in order to be able to complete our comparison to Mindlin and Tiersten’s boundary conditions, is to prove that 


[anti[m • n\ • = [anti[(l — n 0 n) • rh • n] • z/] 

holds also true pointwise. To this aim, let us first remark that 

[anti[(n0 n) ■ rh ■ n]]^ • = [anti[n(n,m • n)]]^ = (n, [rh ■ n] =t: )[anti(n)] =l= • 

= (n, [rh • n]±)(n x = (n, [rh • n]^) t 1 ^, 

where = n x is tangent to curves Or, according to the orientation on 0f2 \ T and T, respectively. 

Hence, we deduce 

[anti[(n 0 n) • rh • n]] + • + [anti[(n 0 n) • rh ■ n]] — • v~ = (n, [rh • n] + ) r + + (n, [rh ■ n\~) r~ 

Therefore, we have 

[[anti[m • n] • v\ = ([anti[m • n]]”*~ — [anti[ra • n]] — ) • v 

= ([anti[(l — n 0 n) • rh • n]]”*~ — [anti[(l — n 0 n) • rh • n]] _ ) • v 
— ((n, [rh ■ n]~*~) r + + (n, [rh • n ]~) r — ) 

= [anti[(l — n 0 n)m • n] • i/] — ((n, [m • n]”*~) r + + (n, [m • n] _ ) r - ). 
Using Stokes theorem and the divergence theorem, we obtain 


/ ((n,[rn • n] + ) , 5u) ds + / ((n, [m ■ n] )r ,Su)ds 

dar dar 

= / (r + , (n, [m • n]”*") <5n) ds + / (r — , (n, [m • n] _ ) du) ds 

dar dar 

= / (n, curl ((n, m • n) du)) da + / (n, curl ((n, rh ■ n) 5u)) da 

JdQ\ r dr 

= / (n, curl ((n, m • n) dw)) da = / div[curl ((n, rh ■ n) <5w)] dv = 0. 

JdQ Jdn 


Therefore, we have 


for all 8u E We choose 


and we obtain 


/ ([anti[m • n\ ■ v\ — [anti[(l — n 0 n) • m ■ n] • v\ , dw) ds = 0, 
dar 

8u = [anti[m • n] • v\ — [anti[(l — n 0 n) • rh • n] • i/J 

/ || [anti[m • n] • i/] — [anti[(l — n 0 n) • m • n] • i /]|| 2 ds = 0. 

dar 

Let us consider an arbitrary parametrization 7 : [a, 6] —y dT of the curve Or. We obtain 

J ^||[anti[m • n] • v\ — [anti[(l — n 0 n) • rh • n\ • ^]||(7(s)) |7 7 (s)|^ ds = 0, 


which implies 


Hence 


which is equivalent to 


||[anti[m • n] • uj — [anti[(l — n 0 n) ■ rh ■ n] • ^]||( 7 (s)) = 0, Vs E [a, 6]. 
|| [anti[m • n] • z/] — [anti[(l — n 0 n) • rh • n] • 1 /] || = 0 on Or 
[anti[m • n] • v\ = [anti[(l — n 0 n) • m • n] • 1 /] on Or. 


(A.32) 


(A.33) 


(A.34) 


(A.35) 


(A.36) 


(A.37) 

(A.38) 
(A.39) 


(A.40) 
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A.7 The proof of Proposition |4.1| 

First, we prove the following lemma: 

Lemma A.l. Let us consider a smooth level surface E = {(xi i X 2 ,xs) E R 3 | F(x\, # 2 , x%) = 0}, 

F : uj C R 3 —> R 3 of class C 2 , then it holds true: 

V[anti(n)] :T = 0. (A.41) 


Proof. Let us first remark that V[anti(n)] : T is well defined, since only the tangential derivative of anti(n) are involved. The 
normal vector to the surface E is given by 


VF 

“ liVFjp 


(ni,n 2 ,n 3 ) 


1 


(Fl,F 2 ,-F 3 ). 


(A.42) 


Let us remark that 


[V[anti(n)] : T] . 


[anti(n)]ij ifc T jk = 



Tjk = [[|VF|| 1 anti(Vi 7 ’)]ij i fc Tj k 

ij,k 


[IV^’ll _1 ],fc[ a nti(VF)]ij- Tj k + HVFir 1 [anti(VF)] iiifc T jk . 


(A.43) 


We compute 

[IIVFir 1 ],* = [(FiF,i)" 1/2 ],fc = [(F,iFi)- 3 / 2 ]2(FiF lk ) = -||VF||- 3 Fj F lk , (A.44) 

anti {n) ijk = -[UjiF^]^ = -eijiFj k , 

Tjk = Sjk - rijn k = Sjk - \\'VF\\~ 2 FjF^ ■ 

Moreover, using ( |A. 10| |, we have 

[anti(VF)]y T jk = [anti(VF) ■ T\ ik = ||VF|| [anti(n) ■ T] ik 

= ||VF|| [anti(n) • (1 - n® n)] ik = ||VF|| [anti(n)] jfc = -||VF|| e ikl rn = -e ik iFj. (A.45) 

Using | |A.44} | |A.45[ in ( |A.43| ) , we obtain 

[V[anti(n)l : T]. = (-||VF||- 3 F,i F, lk ) (-e iks F, s ) + IIVFII” 1 {-e ijt Fj k ) (S jk - || VF||- 2 F,-F fc ) 

= IV-F’ll -3 ( e iks F,s Fj Fj k + eijiF t i k FjF' k ) - || VF|| _1 €ijiF^ k 5j k (A. 46) 

= l|VF|r 3 {e ik j Fj F^Fjk + £ijiF,i k FjF k ) - || VF|| _1 e ik iFj k 
= l|VF||- 3 (euj F,j F ik F k i + eijiF iik FjF ik ) 

= ||VF|| -3 (euj + £iji) FjkFjF'k = 0, 

and the proof is complete. □ 


We now proce ed w ith the proof of Pr oposit ion \^.1\ We start remarking that, using the properties of the projectors T and Q 
introduced in (M2} and the definition ( |3.39| ) of the tensor B one has 

B = -anti (rh • n) = —anti ((T + Q) ■ rh ■ n) = -anti (T • rh • n) + —anti (Q • rh • n) 

= - anti (T • rh ■ n) + — anti (n® n • rh ■ n) 

= -anti (T • rh • n) + —anti ((n, (sym rh) ■ n)n) =: A + —anti (-0 n ). 

Recalling the definition EH of the anti-operator, the second term in the traction boundary condition ( |3.43| ) can hence be manip¬ 
ulated as follows: 


• T) : T]^ — yAijTjp ^ijk' i P n k^jp^j ^hp — (AijTjp) ^Tftp 2 [fijk'4 ,ri k r ^jp) ^Thp- 
We can now remark that 

(fijk'4* ft Tfip = €ijk \fP ,h,Tjh n k H - ^^jh^k^h -^jp,h.2~ftp7lfc] • 

On the other hand, recalling that T + Q = 1 and that Q = n <g> n is a symmetric tensor, we also have 


0 — (Tjp + Qjp\h — Fjp^h + T Up^Uj =y- Tjp,h — 27?.^ft,Tip. 

Using this last equality in ( |A.48| ) and the fact that T ■ n = 0 we finally have 

{fijk^P n kTjp') ^ Tftp = €ijk jh^jh^k T ^ r Tjh n k,h 2"0 71j,ft77.pTftp?7.ftJ 


— ^ijk'^P ihTjh'W'k 4 “ €ijk'^ r Tjh n k,h 0 ? 


(A.47) 
(A.48) 
(A.49) 


(A.50) 
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and we can hence rewrite the term in | |A.47| as 

[V (B • T ) : T] i = yAijTj p — Thp 

= (-A-ijTjp)^ Thp ~ (tijk'lft ,hTjh n k "t" ^ijk^ Tjh n k,h) 

= (AijTjp) h Thp + ~ (tijk'lft ,h(3jh Tjh $jh) n k) ~ (^ijk'fi '^'jh l ^'k,h) 

= (AijTjp’)^ T^p + — {(■ijk'^P ,h' r ^j n h n k) ~ (tijk'tft ,j n k) ~ (fijkrf h^k,h) 
= {AijTjp) h Thp + 0 — {tijk'lft ,j n k) — {tijk'lfr r Tjh n k,h) 

= {A-ijTjp) ^ Thp ~ {^-ijk^ ,j n k) ~ {^ijk^P ^jh^k^h) 
or equivalently in compact form (see the definition of vector product in 1Hh 

V (B -T) :T = V (A-T) :T - i(V^) X n - ^if> V[anti(ra)] : T 

= V (A ■ T) : T + A X (Vij> • T) - V[anti(n)] : T = 

= — V [anti (T • fh • n) • T] : T + —n x [V ((n, (symm) • n))] — — 0 V[anti(ra)] : T. 


(A.51) 


Moreover, using Lemma |A.l[ we have 
Therefore, the proof is complete. 


V[anti(n)] :T = 0. 


A.8 Some lemmas useful to understand Mindlin and Tiersten’s approach 

Lemma A.2. Let Q be an open subset o/M 3 , T an open subsets of dtt and t,g : dLl \ F —> R 3 two functions. Then, the 

y (t, 5u) da + / (g, (1 — n <g> n) ■ curl Su) da = 0 

dn\T JdQ \r 

holds for all Su E C 2 (r2) if and only ift = 0 and (1 — n (g) n) _ = 0. 

Proof. Similar calculations as in Appendix lead to 

/ (£, Su) da + / (< 7 , (1 — n (g) n) • curl da 

JdCi\T JdCi\T 

= (t — V[(anti[(l — n <g> n) - 7/]) • T] : T, Su) da 

Jd n\T 

— / ((1 — n <S> n) • anti[(l — n® n) • 5 ] • n, V<5w • n) da 

JdQ\T 

— / ([anti[(l — n (g) n) • g] • v\, Su) ds 

JddQ\r 

for all variations Su E C 2 (Q). Therefore 

y (t, Su) da + (g, (1 — n <S> n) • curl Su) da = 0 

dQ\T JdQ\T 

implies 

f ( t -V[(anti[(l — n (g> n) • 5]) • T] : T, <5u) da 

Jdn\T 2 

— / ((1 — n <S> n) • anti[(l — n (g) n) • <j] • n, \7Su • n) da 

JdQ \r 

— / ([anti[(l — n 0 n) • 5 ] • i/], <5u) ds = 0 

</ddn\r 


(A.52) 


(A.53) 


equality 
(A.54) 


(A.55) 


(A.56) 


lddci\r 

for all variations Su E C 2 (Q). Since E C 2 (ft), this fact is equivalent to 


(A.57) 


/ (t — V[(anti[(l — n (g> n) • ^]) • T] : T, 5u) da 
Jdd\T 

— / ((1 — n 0 n) • anti[(l — n <S> n) • g\ • n, \7Su • n) da = 0 

Jdn\T ^---—v--- / 


(A.58) 


completely <5u-independent second 
order normal variation of gradient 
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for all variations Su E C 2 (£l). Having now the possibility to consider independent variations of Su and \75u-n on <9f2\T, we obtain 
that 


t — V[(anti[(l — n 0 n) • ol) • T] : T _ = 0, 

lan\r 

(1 — n (g> n) ■ antilYl — n <g> n) • ol • n _ = 0 . 

\dQ\r 

We also remark that 

(1 — n 0 n) • anti[(l — n (g) n) • g\ ■ n = anti[(l — n (g) n) • 'g] • n. 

Hence, it remains to prove that from 

t — V[anti[(l — n (g) ft) • gl • T] : t\ =0, 

L LV J J la^\r 

anti[(l — n ® n) • o| • n _ = 0 

lan\r 

it follows that t\ _ = 0 and (1 — n (g) n) • q\ _ = 0 . 

\dQ\r Iaf2\r 

In the suitable local coordinate system (r, v, n) at the boundary we obtain 

n (g) n = diag( 0 , 0 , 1 ), (1 — n (g) n) = diag(l, 1 , 0 ), g = g T r + g u v + g n n , 

(1 — n (g) n)g = g T r -\-g u u, [(1 — n® n) -g\ x n = g T v — g u t. 

Hence 

anti[(l — n <g) ft) • < 7 ] • n _ = ( g T u — ^ r) _ = 0 


(A.59) 


(A.60) 


(A.61) 


ian\r 


aci \r 

g T _ = 0, _ = 0 (1 — n ® ft) • <7 _ = 0. 


\dQ\r 


laQ\r 


an\r 


(A.62) 


(A.63) 


Moreover, we deduce 


{V[anti[(l — n (g) n) • g\ • T] : T}J _ = [anti[(l - n <g) n) • < 7 ]^ T jk \ _ = [anti[(l - n (g) n) • 5]^-^ r fc 7 y 


la^\r 


a^\r 


lan\r 


la^\r 


= (V{[anti[(l — n (g) n) • < 7 ]^}, r) 

since [anti[(l — n 0 n) • gl« 7 - _ = 0, for all i, j = 1, 2, 3. Therefore, it follows that also t\ _ = 0. 

J J \dQ\r lamr 


the equality 


= o, 

r 

(A.64) 

. *1 =0. 

□ 

Isn\r 


R 3 , g : dfl\T 

—> R 3 two functions. Then, 

= 0 

(A.65) 


holds for all Su E C 2 (Q), <5u|^ = 0 and (1 — n (g) n) • curlSu 


/ (t, <5u) da + / (< 7 , (1 — n (g) n) • curl <5u) da = 0 

ia^\r v/r 

_ = 0 if and only iftl _ = 0 and (1 — n<g> n)q\ = 0 . 

an\r lan\r lr 


Proof. The proof is similar with the proof of Lemma | A.2 


□ 


A.9 Concluding diagrams 
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Standard boundary conditions in the indeterminate couple stress model II59II 

Geometric (essential) boundary conditions (3+2) [only weakly independent] 


(1 — n<S> n) ■ Vu • nL — (1 — n <g> n) ■ 


or (1 — n <S> n) • curl u\ r = (1 — n (g) n) • 6 ext 
t = Divm = | anti(Divra) G so(3) 


Mechanical (traction) boundary conditions (3+2) 

((ct - r) ■ n - \n x V[(n, (symm) • n)]) | 0fAr = t ext , 

(1 — n 0 n) • rh ■ = (1 — n (g> n) • m ext 

Boundary virtual work 

— f dQ ((a — t) • n,8u) da — • n, axl(skew \78u)) da = 0 

— — r) • n — x V[(n, (symm) • n)] J, 8u ) da — • n, |(1 — n (g> n) • [axl(skew V<hi)]j-) da = 0 


3 be 
2 be 


normal curvature 


$ 


Standard boundary conditions in the indeterminate couple stress model, index-format 

Geometric (essential) boundary conditions (3+2) [only weakly independent] 

+ r = “T* e R3 ’ ( £ikm,k - e jkl u l,k n j n i) | r = <HklUl*k - e jkl u ^Jk n j n i! 

or (ui, k n k - u jtk n k 7ijni ) | r = uf *£n k - u e ^n k rijni 


_ _ /.ext 

\dci\r ~ L i ’ 

0 Un\r = m t xt - 


7 ij — o ^ S0(3) 


3 be 
2 be 


Mechanical (traction) boundary conditions (3+2) 

(jTlik n k m jk n k' n 'j n i') 

Boundary virtual work 

— Jqq ({cij — Tij) rij) 8ui da — f@Q(fn • n, axlfskew V<5u)) da = 0 

Sdfl ) +7 2 ^ikl'^k(yO'ij'ft'i'^j > ),l) da ^ Jqq (j^ik^k '^’jk^'k^'j^'i) {^■ikl^'^'l,k ^■jkl^'^ , l,k'^ , j'^ / i) da = 0 


Figure 2: The standard boundary conditions in the indeterminate couple stress model which have been employed 
hitherto by all authors to our knowledge. The virtual displacement is denoted by Su G C°°(r2). The number of 
traction boundary conditions is correct, but the split into independent variations at the boundary is no taken 
to it’s logical end. 
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Boundary conditions in the indeterminate couple stress model in terms of gradient elasticity 


and third order moment tensors 

Geometric (essential) boundary conditions (3+2) [strongly independent] 

w| r = u ext E R 3 , (1 — n 0 n) • Vu • n| r = (1 — n 0 n) • a ext E R 3 , or (1 — n (g> n) • curl u\ T = (1 — n <g> n) • 6 ext 

Mechanical (traction) boundary conditions (3+2) 

((cr — Divm) • n — V[(m • n) • (1 — n ® n)] : (1 — n <g> n)) = £ ext , m = Dvvu[^curv(V[axl(skew Vu)])] 3 be 


(1 — n (g) n) • [m • n] • n| af2 yp = (l-n®n)'i 


[(m-n)- i/]| ar = ir ext , 


“edge line force” on dr 


Boundary virtual work 

— Jqq ((cr — Div m) • n, Su) da — (m • n, \7Su) da = 0 <^> 

— /qq((^ — Divm) • n — V[(m • n) • (1 — n <g> n)] : (1 — n <g> n), Su) da — ((1 — n <g> n) • [m • n] • n, \7Su • n)da 

— Jar (IK™' n )' Su ) ds = 0 

J L equivalent 


Boundary conditions (3+2) in the indeterminate couple stress model in terms of gradient elasticity, 
third order moment tensors, and written in indices 

Geometric (essential) boundary conditions (3+2) [strongly independent] 

Ui |p = «f xt , (ui,k n k ~ Uj' k rijn k ni )| r = uf£n k - u^nirijn k , 

or (e M ui tk - ej k iui tk njni)\ t = (e iki uf^ - e jk iuf£rijni), 

Mechanical (traction) boundary conditions (3+2) 

(+7 tn ijk^k') n j (t &ipk n k rrl ijk hl k Uj 'l.p ) ^ (fifth hip, lip j j fjC'p [' = ^f Xt ) *Tl [j k = D u ^ p 11+u rv ( (+/ k ' ll k.} ) m ) ' 


nijpTLj m p jk n k njrii) hip | a ^ |. — m 

[m p j k n k y,]| ar = + 


_ r! p _ mr *n„n 


p xt n p ni 


“edge line force” on i)T 


Figure 3: The standard strongly independent boundary conditions in the indeterminate couple stress model 
in terms of a third order couple stress tensor coming from full gradient elasticity. The virtual displacement is 
denoted by Su £ C°°(Q). The summation convention was used in index notations. 
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Strongly independent boundary conditions in the indeterminate couple stress model ^ 

Geometric (essential) boundary conditions (3+2) 

n| r = u ext G M 3 , (1 — n 0 n) • Vu • n| r = (1 — n 0 n) • a ext G M 3 , or (1 — n 0 n) • curlu| r = (1 — n 0 n) • 6 ext 
Mechanical (traction) boundary conditions (3+2) 

^(cr — r) • n — 1 n x V[(n, (sym rh) • n)] 

— 1 V[(anti[(l — n 0 n) • rh • n]) • (1 — n 0 n)] : (1-n® n)^ = i ext , 3 be 

(1 — n 0 n) • anti[(l — n 0 n) • rh • n\ • ^I^Qyp — (1 — n 0 n) • m ext 2 be 

[anti[m * ri] * | = 7r ext , “edge line force” on dr 3 be 

Boundary virtual work 

— /qq((o' — t) - n, Su) da — fg^(hh • n, axl(skew X75u)) da = 0 O 

— J^(| (o — t) • n — ^ n x V[(n, (sym rh) • n)] Su) da + fg^ihS • n, |(1 — n 0 n) • [axl(skew V<hx)]|) da = 0 

V 

normal curvature 

— fg^((o — t). n — ^n x V[(n, (sym rh) • n)] — ^V[(anti[(l — n<g> n) • rh • n]) • (1 — n 0 n)} : (1 — n 0 n), <hx) da 

~ \ ~ n ) ' anti[(l — n 0 n) ■ rh • n] ■ n, X7Su ■ n)da — ^ /^ r ([anti[m • n] • i/], Su)ds = 0 J 




equivalent 


Equivalent strongly independent boundary conditions in the indeterminate couple stress model 

Geometric (essential) boundary conditions (3+2) 

n| r = w ext G M 3 , (1 — n 0 n) • Vit • n| r = (1 — n 0 n) • a ext G M 3 , or (1 — n 0 n) • curlu| r = (1 — n 0 n) • 6 ext 
Mechanical (traction) boundary conditions (3+2) 

((cr -r)-n — | V[(anti(m • n)) • (1 — n 0 n)] : (1 — n 0 n)) = t ext , 

(1 — n 0 n) • anti[m • n] • = (1 — n 0 n) • m ext 

[anti[m • n\ • ^]|^ r = 7r ext , “edge line force” on dr 

Boundary virtual work 

— fgfi((o — t) • n , Su) da — fg^(rh ■ n , axl(skew V<5ri)) da = 0 <^> 

— fgft((o — t) • n — ^V[(anti(ra • n)) • (1 — n 0 n)] : (1 — n 0 n), <5it) da 

~ \ — n ® n ) anti(m ■ n) • n, \7Su ■ n)da — ^ /^p([anti[m • n\ • i/J, Su)ds = 0 




equivalent 


3 be 

2 be 

3 be 


Alternative equivalent strongly independent boundary conditions, index-format 

Geometric (essential) boundary conditions (3+2) 

u + = “! xt e r3 > {Ukm,k - tjkiui^njrii) | r = a| xt 

or (ui t knk - Uj }k nkiijni) | r = 6? xt 
Mechanical (traction) boundary conditions (3+2) 

((<Tij -Tij)nj + \(e ihk rh ks n 3 - eij k rhk s n s njn h ), P (Sh P - n h n p )) | an \ r = i| xt , 3 be 

( e ipk m ka ns - ej pk fh ks n a njm) n P | an \ r = m f Xt _ m^n p ni, 2 be 

l^ipkf^ks^s ^pllgp = ’ “edge line force” on <)T 3 be 


Figure 4: The possible boundary conditions in the indeterminate couple stress model. The equivalence of the 
geometric boundary condition is clear. The virtual displacement is denoted by 5u € C°°(Q). 
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